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EXECUTIVE  SUMMARY 


This  report  covers  work  accomplished  at  the  Optical  Sciences  Center  during  the  period  from 
October  1992  through  November  1993,  and  supported  by  the  Joint  Services  Optics  Program.  The  report 
describes  results  obtained  in  1 1  research  projects  spanning  the  field  of  optics  from  fundamental  quantum 
optics  theory  to  applied  optical  engineering.  These  projects  supported  more  than  24  graduate  students 
and  several  post-doctoral  associates,  and  resulted  in  more  than  127  publications  during  this  year.  The 
research  involved  collaborations  at  a  number  of  Department  of  Defense  and  Department  of  Energy 
laboratories  including  the  Army  Center  for  Night  Vision  and  Electro-Optics,  Rome  Air  Development 
Center,  Air  Force  Phillips  Laboratory,  Lawrence  Livermore  National  Laboratory,  and  Sandia  National 
Laboratory. 

Some  of  the  highlights  of  the  research  are  listed  below  and  the  details  of  the  work  are  given  in 
the  following  sections. 

*  A  simulation  model  was  developed  for  a  multiple  quantum  well  infrared  photodetector  that 
enhances  the  ability  to  design  and  optimize  new  detector  structures. 

*  An  improved  near-field  scanning  optical  microscope  was  developed  and  tested. 

*  High  quality  single-crystal  layers  of  beryllium  were  grown  on  germanium  by  molecular  beam 
epitaxy  (MBE)  techniques. 

*  The  properties  of  vertical-cavity  surface-emitting  lasers  (VCSELs)  with  injection  seeding  were 
characterized. 

*  Photon  correlation  measurements  were  used  to  show  that  there  is  no  thresholdless  laser. 

*  Many-body  effects  for  two-component  plasmas  were  characterized  for  quantum  well  structures. 

*  Femtosecond  spectral  hole-burning  techniques  were  employed  to  characterize  the  relaxation  of 
nonequilibrium  carriers  in  semiconductors. 

*  The  properties  of  microcavity  lasers  were  analyzed  and  found  to  be  significantly  different  than 
normal  lasers. 

*  Scanning  tunneling  microscope  techniques  were  developed  to  characterize  semiconductor  optics 
on  a  nanosecond  time  scale  with  nanometer  resolution. 

*  A  figure  of  merit  was  developed  for  the  use  of  quantum  dots  as  all  optical  switches  and  applied 
to  a  variety  of  materials. 

*  Theoretical  models  were  developed  to  characterize  quantized  atomic  motion  in  one-,  two-,  and 
three-dimensional  optical  lattices. 

Richard  C.  Powell,  Director 
Optical  Sciences  Center 
February  1994 
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RESEARCH  FINDINGS 

A  proposal  for  a  simulation  model  of  a  multiple  quantum  well  infrared  photodetector  (MQWIP) 
is  presented.  Although  the  theory  of  MQWIPs  is  well  studied,  very  little  simulation  results  of  MQWIP 
have  been  published.  Besides  understanding  the  operaion  of  MQWIPs,  the  simulation  nKxlel  is  a  useful 
tool  for  designing  and  optimizing  new  detector  structures. 

The  theoretical  analysis  is  based  on  Boltzmann  transport  equation,  which  is  solved  with  Monte 
Carlo  (MC)  simulation  combined  with  a  Poisson  equation  solver.  The  MC  method  applied  to  charge 
transport  in  deteaor  consists  of  the  simulation  of  carriers  in  a  given  electric  field  and  scattering 
mechanisms.  The  approach  is  semiclassical  electron  trajectories  and  quantum  mechanically  calculated 
scattering  probabilities. 

A  widely  used  method  of  modeling  semiconductor  devices  is  the  drift  diffusion  (DD)  model, 
which  uses  phenomenological  parameters  such  as  mobility  and  diffusion  coefficient.  It  assumes  that  the 
energy  and  momentum  of  the  charger  carriers  remain  close  to  that  of  an  equilibrium  distribution.  This 
is  a  good  approximation  for  many  large-scale  devices;  however,  in  submicron  range  this  assumption  is 
no  longer  v^id.  For  multiple  quantum  well  (MQW)  structures  a  more  elaborate  method  is  asked  for. 
MC  methods  give  self-consistent  solutions  for  the  Boltzmann  transport  equation  which  gives  a  better 
description  of  multiple  quantum  well  devices. 

Monte  Carlo  Simulation 

The  MC  method  applied  to  charge  transport  in  semiconductors  consists  of  simulating  one  or  more 
carriers  in  a  given  electric  field  and  scattering  mechanism.  The  simulation  starts  with  approximating  the 
electron  charge  density  This  could  be  calculated,  for  example,  with  the  DD-model.  This  is  used  as 
input  to  the  Poisson  equation. 

V*E  -  p„  (1) 

where  p^  is  the  total  charge  density  from  electrons  and  ionized  donors.  The  MC  part  consists  of 
simulated  electron  flights:  one  electron  is  given  initial  conditions  with  wavevector  k.  Then  the  duration 
of  the  first  free  flight  is  chosen  with  a  probability  distribution  d^ermined  from  the  scattering 
probabilities.  By  keeping  statistics  of  the  electron  motion  (i.e.,  how  much  time  it  spends  in  each  lattice 
point),  it  is  possible  to  estimate  the  charge  density  inside  the  detector.  After  many  simulated  trajectories, 
the  average  electron  charge  density  p,,  is  calculated  and  fed  back  to  solve  the  Poisson  Eq.  (1)  again.  The 
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entire  process  is  repeated  until  a  self-consistency  between  the  electric  field  and  the  charge  density  is 
acquir^. 

During  the  free  flight  the  electric  fields  are  acting  according  to  the  relation 

Ak  -  -eE  (2) 


and  kq>t  constant.  All  the  quantities  of  interest,  velocity,  energy,  etc.,  are  recorded.  A  scattering 
mechanism  is  chosen  at  the  end  of  the  free  flight.  New  state  of  k  is  calculated  according  to  the  chosen 
scattering  mechanism  and  it  becomes  the  initial  state  for  the  new  free  flight. 

Since  the  electron  scattering  probabilities  are  functions  of  electron  energies,  calculation  of  the 
flight  durations  would  be  difficult  and  time  consuming.  To  overcome  this,  the  self-scattering  scheme  is 
introduced.  The  maximal  scattering  rate  F  is  used  as  a  constant  to  produce  flight  durations.  After  the 
flight,  it  is  decided  which  physical  scattering  mechanism  to  use  or  whether  to  scatter  at  all.  In  the  latter 
case,  it  is  said  that  the  electron  underwent  a  self  scattering.  If  the  scattering  rates  vary  over  space  and 
k-values  a  lot.  most  of  the  simulated  scattering  events  are  self-scattering.  This  slows  the  simulation 
considerably  and  extensive  algorithm  development  has  been  done  to  And  efficient  ways  to  reduce  the 
amount  of  self  scattering. 

The  result  of  the  calculation  becomes  more  precise  as  the  simulation  goes  on  and  the  simulation 
ends  when  the  quantities  of  interest  are  known  with  the  desired  precision. 

The  MC  method  simulates  the  dynamics  of  a  single  electron  inside  the  crystal,  which  makes  it 
similar  to  experimental  techniques.  In  fact,  it  can  be  considered  as  a  simulated  experiment. 

Scattering  Rates 

For  the  two-dimensional  (2D)  systems  the  scattering  rates  are  calculated  as  for  the  bulk,  for  lack 
of  a  better  model.  Scattering  rate  calculations  usually  rely  on  flrst-order  perturbation  theory.  Fermi's 
golden  rule  is  a  constant  rate  of  transition.  Scattering  in  three  dimensions  (3D)  is  considered  as  well- 
understood  theory,  and  comprehensive  presentations  can  be  found  in  the  literature. 

Usually  Ae  following  scattering  mechanisms  are  included  in  semiconductor  device  modeling; 
acoustic  phonons,  optical  phonons,  ionized  in^urities,  interface  scattering,  and  carrier-carrier  scattering. 
In  high-temperature  region  (even  in  room  tenqierature)  optical  phonon  scattering  is  the  most  important. 
In  low  temperatures  the  ionized  impurity  scattering  dominates. 


Absorption 

The  photon  absorption  in  2D  was  solved  already  in  1957  by  Elliott: 


2y<«>|<i^|^  f 2m‘ 


4E, 

—  8 

fl«-E"+E  B 

1"-^1 

R 

W  J 

^  ‘'coshA 

(3) 


where  E,  is  the  three-dimensional  exciton  Rydberg  given  by  E,  =  e*m*/2c,V.  The  first  term  in  the 
bracket  corresponds  to  the  excitonic  Imes.  The  excitonic  absorption  is  of  great  importance  in  2D  systems 
because  the  binding  energies  are  four  times  stronger  than  in  bulk.  However,  excitons  are  not  of  interest 
in  detector  operation  because  they  do  not  produce  photoelectrons. 

The  absorption  constant  cr^  can  he  calculated  by  solving  the  transition  probability  matrix  d,,, 
where  c(v)  refers  to  conduction  (valence)  band. 
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The  from-first  principles  calculation  of  the  probability  matrix  (1„  is  rather  complicated  and  needs 
super-computer  power.  However,  since  the  steady-state  MC  calculation  already  predicts  the  charge 
density,  we  can  use  this  to  determine  the  absorption  (up  to  a  constant),  because  it  is  proportional  to  the 
charge  distribution.  The  constant  can  be  fixed  using  experimental  data  of  the  absorption  rate. 

The  absorption  in  2D  structures  depends  on  polarization.  (Dnly  the  light  polarized  perpendicular 
to  the  layers  is  absorbed. 

Simulation 

The  simulation  of  MQWIP  can  be  divided  in  two  parts;  steady-state  (or  carrier-transportation) 
and  absorption  simulation.  Figures  1  and  2  are  corresponding  flow  charts.  The  purpose  of  the  steady- 
state  simulation  is  to  calculate  how  the  electrons  and  doped  ions  are  distributed  in  the  medium  (elecuic 
field),  how  fast  the  electrons  reach  the  collector  (mean  velocity)  and  how  many  electrons  reach  the 
collector  (photo-conductive  gain).  From  the  charge  distribution  the  absorption  coefficient  a(z),  where 
z  is  the  spatial  coordinate,  can  be  determined  up  to  a  constant. 

In  absorption  simulation  the  place  z  of  absorption  is  deterinined.  The  initial  energy  and 
polarization  is  given  to  the  photon  entering  the  detector.  Probabilistic  criteria  are  used  to  determine  if 
the  photon  is  reflected  from  the  surfa^  If  not,  a  stmistical  distribution  calculated  from  the  absorption 
equation  e**^  is  used  to  determine  if  the  photon  is  absorbed  to  the  first  quantum  well.  The  same  is 
repeated  for  the  following  wells  until  the  photon  is  absorbed  or  until  it  reaches  the  collector.  Combining 
the  absorption  rate  with  the  simulated  transport  properties  of  the  detector,  the  pulse  shape  and  quantum 
efficiency  can  be  determined. 
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Figure  2.  Flowchut  of  ifae  ebsofpdoa  simulelioa. 
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Data  Storage,  Maui.  H  aii,  July  1993. 

J.L.  Kann,  T.D.  Milster,  F.  Froehlich,  R.W.  Ziolkowski  and  J.  Judkins,  ‘Numerical  analysis  of  a  two 
dimensional  near-field  prtfoe,”  poster  presentation.  Second  Imernational  Conference  on  Near  Field  Optics, 
Raleigh,  North  Carolina,  October  1993. 

J.L.  Kann,  T.D.  Milster,  F.  Froehlich,  R.W.  Ziolkowski  and  J.  Judkins,  'Numerical  analysis  of  a  two 
dimensional  near-field  probe,”  accepted  for  publication,  J.  Ultramicroscopy.  June  1994. 


RESEARCH  HNDINGS 
Experimental  Work 

We  are  continuing  to  develop  our  second-generation  near-field  scanning  optical  microscopy 
(NSOM)  instrument  that  we  reported  on  previously.  The  current  configuration  of  the  instrument  is  shown 
in  Fig.  I,  and  is  similar  to  die  design  desaibed  in  detail  in  last  year’s  rqiott.  The  figure  shows  the  head 
assembly  of  die  microscope  consisting  of  optical  and  mechanical  hardware  with  the  servo  system.  We 
have  interfaced  our  head  assembly  to  a  commercial  tuomic  force  microscope  (AFM)  (not  shown)  that  was 
donated  to  the  project  by  WYKO  Inc.,  Tucson,  Arizona.  The  AFM  is  us^  to  control  the  raster  scanning 
of  the  sample,  and  to  acquire,  di^Iay  and  process  the  dau.  Our  microsct^  unfortunately  is  still  not 
fully  operatiooal  due  to  difficulties  encountered  with  the  servo  system  diat  maintains  the  probe  tip-to- 
sample  spacing  (10  nm  to  40  nm)  required  for  high  spatial  resolution.  The  servo  system  represents  a 
significant  engineering  challenge,  and  we  have  solved  numerous  hardware  problems  during  the  past  year 
that  were  hindering  its  operatioo.  However,  the  sample  stage  continues  to  be  a  problem  due  to  its  poor 
frequency  reqionse  and  resonance  behavior.  This  stage  is  a  complex  piezoelectrically  driven  flexure 
design  diat  provides  raster  .scanning  of  the  sample  in  x  and  y,  and  servo  correction  m  z.  It  is  our  opinion 
that  die  electro-mechanical  performance  of  our  servo  cannot  be  further  improved  to  an  adequate  level 
using  diis  stage  and  our  current  hardware  configuration.  Therefore,  we  are  planning  to  modify  the  head 
of  our  WYKO  AFM  for  use  as  a  NSOM  microscope.  The  WYKO  unit  is  a  proven  design,  and  its  use 
will  eliminate  the  problems  we  are  experiencing  with  our  mechanics  and  piezoelectric  positioners.  The 
cost  of  taking  this  approach  is  sacrificing  the  excellent  scan  linearity  afforded  by  die  flexure  sample  stage, 
and  some  loss  of  system  flexibility.  We  anticipate  carrying  out  the  modificatioos  to  the  WV7CO  unit 
during  the  next  funding  cycle. 

We  are  studying  the  servo  system  in  detail  to  quantify  its  sensitivity  and  to  understand  its 
performance  limitations.  The  servo  system  used  in  our  microscope  is  based  on  sensing  die  idomic  forces 
between  the  tip  and  sample  surfoce  to  maintain  the  tip-to-sample  spacing.  A  small  dither  oscillation  is 
applied  to  the  tip  (typic^ly  several  nanometers),  and  the  amplitude  of  die  dither  decreases  as  the  tip 
comes  close  to  the  surfoce  and  experiences  the  attraaive  van  der  Waals  forces  of  the  surface  atoms.  The 
sensing  of  this  minute  didier  amplitude  is  a  critical  part  of  the  servo  system.  The  sensing  is  performed 
by  diffracting  light  from  a  focuUd  HeNe  laser  off  the  tip  and  by  diffwentially  deteaing  the  resultant 
modulation  of  the  diffraction  pattern  with  a  photodiode.  We  developed  diis  tedmique  simultaneously, 
but  independendy,  with  other  NSOM  research  groups.  This  technique  is  superior  to  the  scheme  we 
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originally  proposed  in  last  year’s  report,  which  involves  imaging  the  dithering  tip  to  a  pinhole.  We  are 
conducting  a  ^eoretical  and  experimental  analysis  of  this  technique  for  publication,  as  no  such  analysis 
has  yet  appeared  in  the  literature.  Figure  2  shows  dieoretical  calculations  of  the  differential  signal  vs 
lateral  tip  offset  with  defocus  as  a  parameter.  The  strength  of  the  signal  is  proportional  to  the  slope  of 
a  given  curve.  Therefore,  the  modd  indicates  that  the  sensitivity  is  enhanced  by  defocusing  the  tip  from 
the  laser  spot,  and  defocusing  makes  the  signal  less  dependent  on  tip  offset. 

Modeling  Work 

A  rigorous  electromagnetic  modd  is  being  used  to  study  both  the  near-zone  (within  100  nm  of 
the  aperture)  and  far-zone  (out  to  hundreds  of  microns)  fidds  associated  with  a  typical  NSOM  structure. 
A  two-dimensional  ffnite-difference-time-domain  (FDTD)  code  is  used  to  find  the  steady-state  time- 
averaged  fidds  in  and  around  a  sub-wavdength  ^)erture.  The  FDTD  code  discretizes  the  problem  in 
both  time  and  space  (creating  a  grid),  and  then  solves  Ampere’s  and  Faraday’s  equations  within  the  grid 
at  each  time  stq>.  llie  code  runs  for  16  optical  cycles,  and  data  is  time  averaged  over  the  last  cycle. 
This  code  also  allows  us  to  study  the  effects  of  placing  various  multi-layer  structures  in  proximity  to  the 
^>erture.  An  angular  spectrum  technique  is  then  applied  to  propagate  the  fidd  components  from  the 
FDTD  code  to  the  fv-zone. 

An  example  output  from  the  modd  is  shown  in  Fig.  3.  The  figure  shows  a  perspective  plot  of 
the  normalized  dectric  ffdd  intensity  raunating  from  a  slit  aperture.  In  the  modd,  a  normally  incident 
TM^iarized  plane  wave  (1  -  SOO  nm)  is  used  to  illuminate  a  1(X)  nm  wide  slit  in  a  metal  screen  of 
infinite  conductivity.  The  metal  screen  is  25  nm  thick,  and  the  medium  surrounding  it  is  air.  The  fidd 
enhancement  near  the  aperture  is  caused  by  the  sharp  comers  of  the  slit  in  the  metal,  where  a  large 
charge  density  accumulates.  The  energy  past  the  aperture  remains  wdl  collimated  for  a  distance  of 
approximatdy  half  the  a|>erture  width.  Past  this  distance,  it  diverges  linearly,  similar  to  what  is  seen  for 
a  larger  aperture.  However,  the  intensity  falls  off  exponentially,  due  to  die  fact  that  most  of  the  energy 
transmitt^  through  the  aperture  is  not  propagating,  but  evanescent. 
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RESEARCH  FINDINGS 

This  year  we  continued  our  study  of  epitaxial  beryllium  growth  using  molecular  beam  epitaxy 
(MBE).  We  were  successful  in  growing  very  high  quality  single-crystal  boyllium  on  germanium  (1 1 1) 
substrates.  Our  epitaxial  Be  films  grown  on  Ge  are  even  higher  quality  than  our  previous  films  grown 
on  silicon  (1 1 1)  a^  sapphire.  The  quality  of  the  initial  Ge  surface  and  the  substrate  temperature  during 
growth  are  the  most  in^rtant  vari^les  affecting  the  structure  of  the  Be  film.  We  also  have  grown 
certain  Be-Ge  hecerostructures. 

We  continued  our  development  of  non-specular  x  ray  scattering  as  a  multilayer  characterization 
tool.  The  non-specular  behavior  can  be  underst^  in  terms  of  current  theories  of  x  ray  scattering  from 
multilayers  with  correlated  interface  rou^ess.  Detailed  analysis  of  the  evolution  of  interface  roughness 
in  multilayers  using  non-specular  scattering  has  been  carried  out  in  collaboration  with  D.  Steams  at 
Lawrence  Livermore  Laboratory.  We  have  used  this  technique  to  determine  the  evolution  of  the 
roughness  power  spectral  density  from  mterface-to-ioierface  in  the  multilayer  film.  We  believe  this  is 
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the  first  time  such  detailed  infbnnatioo  about  the  roughness  evolution  in  multilayers  has  been  measured 
without  a  synchrotron  source. 

Beryllium  Results 

We  have  previously  succeeded  in  growing  single-crystal  hep  beryllium  on  sapphire  (001)  and  Si 
(111)  substrates  using  MBE.  During  this  period,  we  have  succeeded  in  growing  Be  epitaxially  on  Ge 
(1 1 1).  These  films  are  even  higher  quality  than  the  films  we  grew  previously  on  sapphire  and  Si.  The 
higher  quality  is  probably  due  to  the  smaller  lattice  mismatch  between  Ge  and  Be,  1.0  vs.  3. 1  %  for  Be- 
on-Si  and  3.9%  for  Be-on-sapphire. 

BerylUum  on  Silicon.  Table  1  summarizes  the  characterization  results  for  Be  grown  on  (111) 
Si.  These  structural  parameters  are  discussed  m  more  detail  in  our  previous  technical  reports  and  our 
publications.  We  made  additional  samples  at  T=:300  *C  and  T=4(W  "C  to  verify  the  high  quality  at 
300  ”C  and  morphology  change  at  400  The  films  with  the  highest  quality  are  those  dqiosited  at 
T=3(X)  **€.  The  sample  made  at  T-SOO  ‘C  has  better  crystallinity,  but  is  composed  of  isla^  rather 
than  a  continuous  film.  Scanning  electron  microscope  (SEM)  and  atomic  force  microscope  (AFM)  images 
of  films  grown  at  T-400  ”C,  show  die  hexagonal  island  formations  typical  of  these  higher  T  samples. 
These  islands  are  peculiar  in  that  they  are  typically  large  (-1  iim  diameter),  flat-topped,  roughly 
hexagonal  islands  with  terraced  sides  topped  by  a  small  island  in  the  center.  Because  the  morphology 
changes  so  drastically  between  T=3(X)  *C  and  T=400  "C,  there  may  be  a  surface  roughening  transition 
in  this  temperature  region.  Further  studies  of  the  morphology  as  a  function  of  deposition  temperature 
and  annealing  temperature  are  needed  to  understand  the  driving  forces  behind  this  drastic  morphological 
change. 

Table  1.  CharKterizatioo  results  for  epitaxial  beryllium  films  grown  on  Si  (111)  at  several  substrate 
temperatures  (T).  The  effective  thickness  fiom  ion  bcaun  analysis  (da^)  is  a  measure  of  how  much  beryllium 
was  deposited.  Better  crystalline  quality  is  indicated  by  having  an  x-ray  coherence  length  (L)  nearer  to  the 
film  thickness  (d^j^),  and  by  lower  values  of  rocking  curve  width  (Sf)  and  ion  beam  channeling  ratio  (x). 


TCO 

dllA  (A) 

^lAD  (^) 

L(A) 

MFWHM 

(deg.) 

X 

100 

1870±190 

1630±10 

920 

0.38±0.02 

0.73±0.02 

200 

1970±200 

1770±I0 

1050 

0.32  ±0.02 

0.38±0.02 

300 

1700±170 

1530±10 

1320 

0.24±0.02 

0.33±0.02 

400 

1060±110 

islands 

850 

0.26±0.02 

0.41  ±0.02 

500 

1120±80 

islands 

1560 

0.20±0.02 

0.79  ±0.02 

Beryllium  on  Germanium.  Results  from  RHEED  analysis  of  the  Be  fllms  grown  on  Ge  buffer 
layers  on  Si(lll)  indicated  that  the  Be  grows_with  the  Be[0(X)l]  direction  parallel  to  the  Ge[lll] 
direction.  Also,  the  Be  is  oriented  with  Be(l010]  direction  parallel  to  foe  Ge[110]  direction  and  foe 
Be[l  120]  direction  parallel  to  foe  Ge[2111.  In  addition,  a  surface  superstructure,  probtfoly  O^3xi/3)R30‘ , 
is  observed  on  foe  samples  wifo  T  2  300  *C.  X  ray  diffraction  analysis  and  ion  channeling  results 
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indicate  an  increase  in  crystalline  quality  as  T  is  increased.  However,  samples  deposited  at  T  2  300  'C 
are  visibly  rough.  The  effects  on  cryst^line  quality  and  surface  roughness  of  increasing  the  temperature 
T  can  be  explained  through  surface  mobility  considerations.  Increasing  T  probably  results  in  increased 
mobility  of  the  Be  adatoms  allowing  more  ordered  growth.  However,  above  T  =  300  ’C,  the  surface 
mobility  may  become  great  enough  to  favor  formation  of  large  islands.  Increased  defects  and  surface 
roughness  may  result  as  these  islands  grow  together. 

Analysis  of  RHEED  patterns  for  Be  films  grown  on  Ge  buffers  on  Ge(lll)  substrates  indicate 
that  the  Be  grows  with  the  same  orientation  as  when  it  is  deposited  on  Ge  buffers  on  Si(l  1 1)  substrates. 
However,  in  the  case  of  Be  on  Ge(l  1 1),  we  observe  the  surface  superstructure  on  all  sample^  en  for 
those  with  T  <  300  *C.  X  ray  diffraction  analysis  and  ion  beam  channeling  show  the  highest  iline 

quality  is  attained  at  T  «  300  *C.  As  T  increases  above  300  *C  or  decreases  below  300  *C,  alline 

quality  is  decreased.  All  of  the  samples  deposited  were  visibly  smooth.  Again,  the  effect  on  crystalline 
quality  of  increasing  the  temperature  T  may  be  explained  in  terms  of  adatom  surface  mobility,  with  T  = 
300  **C  resulting  in  a  mobility  that  favors  highly  ordered  growth  of  smooth  thin  films.  In  addition,  the 
crystalline  quality  of  these  films  grown  on  Ge(lll)  was  significaDtly  greater  than  Be  films  grown  on 
Ge(lll)  buffer  layers  on  Si(lll).  We  suggest  diat  the  influence  of  the  substrate  material  on  the 
crystalline  quality  of  the  Be  overlayer  is  a  result  of  differences  in  the  quality  of  the  Ge  buffer  layen 
Atomic  force  microscopy  and  ion  beam  channeling  reveal  that  the  Ge  buffer  layers  dqwsited  on  Ge(l  1 1) 
are  of  much  higher  crystalline  quality  than  the  buffer  layers  grown  on  Si(l  1 1).  Since  defects  in  a  buffer 
layer  tend  to  propagate  through  an  overlayer,  the  poorer  quality  Ge  buffers  on  Si  most  likely  lead  to 
poorer  quality  Be  films. 

Gemunium  on  Beryllium.  RHEED  analysis  of  the  Ge  films  grown  on  Be  buffers  on  Si(l  1 1) 
indicate  that  Ge  grows  epitaxially  for  T  =  425  *C,  but  polycrystalline  for  T  s  300  *C.  This  behavior 
also  may  be  explained  in  terms  of  increased  Ge  adatom  mobUity  at  bjgber  temperatures.  For  the  epitaxial 
samples,  the  Ge  grows  with  Ge( III]  |  Be[0001],  Ge[110]  |  Be[1010],  and  Ge[211]  |  Be[1120].  Xray 
diffraction  rocking-curve  analysis  and  ion  beam  channeling  show  that  the  crystalline  quality  of  these  Ge 
films  is  poor  compared  to  very  good  Ge  films,  such  as  Ge  grown  on  Ge(l  1 1). 

RHEED  analysis  of  Ge  films  grown  on  thin  Be  buffers  on  Ge(lll)  indicate  that  Ge  grows 
epitaxially  at  all  temperatures  investigated.  However,  for  films  grown  at  T  ^  400  *C,  faint  transmission 
diffraaion  patterns  are  observed.  These  transmission  patterns  are  characteristic  of  rough  films.  RHEED 
patterns  of  films  grown  at  T  »  500  ’C  indicate  these  films  are  smooth  compared  to  those  grown  at  lower 
temperature  T.  X  ray  diffraction  analysis,  RBS  analysis,  and  ion  beam  channeling  results  are 
inconclusive  for  these  samples  due  to  the  sample  geometry.  Specifically,  the  Be  film  is  too  thin  to  be 
characterized  by  any  of  these  techniques.  Also,  because  the  Be  film  thickness  is  only  -50  A, 
characteristics  of  the  Ge  overlayer  cannot  be  resolved  from  those  of  the  Ge  substrate.  Transmission 
electron  microscopy  CTEM)  is  currently  in  progress  to  characterize  these  samples. 

Results  for  Be  films  grown  on  Ge(l  1 1)  and  Ge  films  grown  on  500  A  ^  buffer  layers  on  Si(1 1 1) 
are  shown  in  Table  2  and  Table  3. 

In  conclusion,  our  results  show  that  very  high  quality  Be  thin  films  may  be  grown  on  Ge(l  1 1). 
We  have  demonstrated  the  effects  of  substrate  temperature  and  substrate  quality  on  the  quality  of  the  Be 
film  and  postulated  some  explanations  for  these  effects.  In  addition,  we  hs  shown  epitaxial  growth  of 
Ge  on  B^OOOl).  This  is  a  significant  first  st^  toward  Be/Ge  heterostructures.  Such  structures  may  be 
potentially  usefiU  in  die  areas  of  x  ray  optics  and  electronics. 
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Table  2.  X-ny  difFn^tioo  and  ioo  beam  channeling  reaultt  for  Be  on  Ce  buffer  layers. 


Substrate 

TCO 

Rocking  Curve 

FWHM 

Xm 

Si(lll) 

200 

.50 

.53 

Si(lll) 

300 

.38 

.52 

Si(lll) 

400 

.35 

.45 

Ge(lll) 

200 

.30 

.25 

Ge(lll) 

300 

.19 

.18 

Ge(lll) 

400 

.20 

.39 

Table  3.  X-ray  difliractioo  and  ion  beam  channeling  results  for  Ge  on  SOO  A  Be  buffer  layers. 


Rocking  Curve 

Substrate 

TCQ 

FWHM 

X» 

Si(lll) 

300 

polyciyiulline 

polyciyftaUine 

Si(lll) 

425 

.60 

.64 

Non-specular  X  ray  Scattering  flrom  Multilayer  X  ray  Mirrors 

We  have  observed  resonant  non-specular  x  ray  scattering  effects  in  Si/Mo  multilayers  and  have 
described  the  experimental  techniques  and  results  in  previous  reports.  Using  a  simple  analysis  of  the 
data,  based  on  recent  theories  of  x  ray  scatt^g  from  periodic  multilayers,  we  previously  used  our 
measurements  to  determine  diat  correlated  roughness  was  present  in  these  multilayers  for  spatial  periods 
/> 200  A.  In  order  to  extract  more  structural  information,  a  specific  model  must  be  chosen.  The  model 
must  mclude  the  power  spectra  of  die  interfaces  and  describe  how  the  roughness  is  propagated  through 
die  stack.  In  collaboration  with  Daniel  Steams  at  Lawrence  Livermore  Laboratory,  fois  type  of  detailed 
analysis  of  die  data  for  a  series  of  Si/Mo  multilayer  (sputter-deposited  at  3,  6,  9,  and  12  mtorr  Ar 
pressure)  has  now  been  completed.  Using  a  growfo  m<xlel  developed  by  Steams  wh«e  die  growth  is 
parameterized  in  terms  of  a  rdaxation  parameter  ¥  describing  the  damping  of  features  widi  increasing  film 
thickness  and  0  representing  the  volume  of  a  growth  unit.  The  results  are  die  power  spectral  density 
(PSD)  of  the  roughness  as  a  function  of  ^atial  frequency  and  number  of  layers  and  a  correlation  number 
as  a  function  of  spatial  frequency.  The  results  show  that  roughness  widi  />  1000  A  is  correlated  over 
hundreds  of  bilayers,  but  die  correlation  number  drops  rapidly  with  increasing  spatial  frequency.  For 
/<  100  A  the  roughness  is  effectively  uncorrelated  (correlation  number  <2  bilayers). 
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Results 

Figure  1  shows  a  measured  scatter  curve  through  the  fifth  order  Bragg  peak  for  a  Si/Mo 
multilayer  with  a  period  of  124.7  A  and  the  calculated  fits.  In  this  figure  the  scattering  angle  is  the  angle 
of  the  detector  with  respect  to  normal.  The  peak  positions  observed  in  this  scatter  curve  correspond  to 
/>200  A.  Calculations  were  performed  with  and  wifoout  die  roughness  of  the  superpolished  substrate 
(ff  » 1  A)  included  in  the  noodel.  In  this  case,  the  substrate  roughness  had  little  effect.  For  samples  made 
on  float  glass  (<r«5  A),  the  substrate  roughness  affected  the  non-specular  intensity  for  peaks 
corresponding  to  s^proximately  />  1000  A. 


Figure  1.  A  measured  curve  through  the  third  order  Bragg  peak  for  a  Si/Mo  multilayer  with  a  period  of 
124.7  A  and  the  calculated  fits.  The  scattering  angle  is  the  positioo  of  the  detector  iiom  normal. 

Good  fits  were  obtained  for  the  3  and  6  mtorr  samples.  Althou^  die  roughness  increased  from 
3  to  6  mtorr,  it  was  at  9  mtorr  diat  the  multilayer  structure  became  severely  distorted  by  columnar-like 
growth.  The  fits  to  the  9  mtorr  are  much  more  uncertain  because  of  the  additional  peak  broadening, 
evident  in  die  sqpecular  data,  diat  is  not  correcdy  accounted  for  by  the  model.  The  results  of  the  fitting 
are  shown  in  Table  4.  The  variation  of  the  parameters  with  Ar  pressure  is  consistent  with  a  simple 
mterpretation.  Recall  that  f  characterizes  the  range  of  adatoms  on  the  surface  and  Q  characterizes  the 
average  cluster  size.  At  higher  pressure,  the  sputtered  atoms  undergo  more  frequent  collisions  on  the  way 
to  the  growing  film,  and  therefore  arrive  with  lower  energy.  Thus,  at  hi^er  pressure,  both  v  and  Q 
might  be  expected  to  decrease  due  to  reduced  surface  mobility,  as  is  seen  in  Table  4. 

In  summary,  we  have  made  non-specular  x  ray  scattering  measurements  on  Si/Mo  multilayers 
using  a  low-angle  diffractometer  and  Cu  K,  radiation.  These  data  have  been  analyzed  quantitatively, 
using  a  dieory  of  x  ray  scattering  from  non-ideal  multilayers  and  a  simple  model  of  thin  film  growth. 
The  analysis  determined  the  (frequency  dependent)  PSD  of  the  interfaces,  and  values  for  parameters  that 
describe  the  film  growth.  The  variation  of  these  parameters  with  pressure  reflects  the  changing  mobility 
of  the  depositing  atoms.  We  now  have  developed  this  technique  to  the  point  that  it  can  be  used  to 
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determine  not  only  the  magnitude  of  the  roughness  but  its  origin.  This  information  is  invaluable  in 
improving  both  the  substrate  surface  roughness  and  the  intrinsic  multilayer  roughness. 


Table  4.  Values  for  the  average  volume  of  the  growth  unit  0.  and  the  relaxation  parameter  r.  as  determined  from 
fits  to  the  acatter  curves. 


Ar  Pressure  (mtorr) 

f(A) 

Q(A^ 

3 

20 

200 

6 

10 

100 

9 

10 

10 

14 
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RESEARCH  FINDINGS 
RoutC’to-Iitjection  Locking* 

We  injected  a  stable  coherent  beam  into  a  vertical-cavity  surface-emitting  laser  (VCSEL)  and 
studied  the  route-to-injection  locking  as  a  function  of  the  injected  power.  Pushing  or  pulling  of  the 
original  lasing  was  seen,  acconqianied  by  the  tfipearance  of  new  frequencies.  The  bdiavior  depends  upon 
the  sign  of  the  detuning  of  the  injection  from  the  original  lasing.  We  have  compared  these  results  with 
a  semiconductor  theory  of  injection  locking  and  have  obtained  excellent  agreement.  This  understanding 
and  ability  to  accurately  model  can  be  very  useful  in  injection-locking  arrays  of  VCSELs. 

By  comparing  the  locking  range  for  positive  and  negative  detuning  A,  it  is  possible  to  estimate 
the  alpha  parameter  a.  Figure  1  shows  the  experimental  data  points  for  reciprocal  locking  power  as  a 
fonaion  of  frequency  detuning,  from  which  we  deduce  an  a  between  1.4  and  2.6.  This  covers  our 
previously  measured  value  2.0  ±  0.5  and  the  plasma-theory  conq)uted  value  of  2.27.* 


Local  Gain  Asymmetry  and  the  Generation  of  New  Frequencies^' 

When  a  stable  coherent  beam  is  injected  into  a  VCSEL  that  is  lasing  just  above  threshold,  we  see 
a  modificatioD  of  the  emission  spectrum  as  the  injected  power  is  increased.  The  emission  is  enhanced 
below  the  injection  frequency  and  reduced  above.  These  few-GHz  modifications  that  occur  for  few- 
microwatt  injections  arise  from  coherent  energy  transfer  involving  spatial  gratings  written  by  the  pump 
and  probe,  which  involve  the  total  carrier  density. 

When  the  injection  is  mid-range  (10  ftW  to  40  /iW)  and  tuned  close  to  the  high-frequency  side 
of  the  VCSEL  lasing  [and  the  VCSEL  lasing  is  strong  and  narrow  1.5  GHz)],  a  new  frequency  is 
generated  at  about  5  GHz  to  the  long  wavelength  side  of  the  injection.  As  foe  injection  power  is 
increased,  foe  new  frequency  shifts  to  longer  wavelengths  by  as  much  as  25  GHz,  lioearly  proportional 
to  the  intracavity  power  at  the  injection  frequency,  see  Fig.  2.  Both  foe  new  frequency  and  foe  injection 
can  have  relaxation  oscillation  sidebands,  showing  that  foe  new  frequency  is  not  a  rdaxation  oscillation 
sideband.  The  distorted  dispersive  shape  of  foe  local  gain  modification  a^  foe  linear  dependence  of  foe 
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Figure  1.  Reciprocal  of  the  injection  power  level  needed  to  frequency  lock  the  VCSEL  as  a  function  of  detuning. 
The  solid  curve  shows  the  functional  form  of  the  locking  range  equation  u  asymmetry  assuming  a  »  2.0. 


InLracavily  Injected  Power  (mW) 

Figure  2.  Acceleration  of  coherent  transfer  of  energy  by  stimulated  emission  and  absorption.  Absolute  detunings, 
\^^\  -  Iw'mI.  firom  two  series  of  spectra  like  that  m  (a),  follow  calculations  of  the  gain  peak  (solid) 
and  dip  (dashed);  from  Ref.  2 
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new  gain  peak  upon  tntracavity  injected  power  are  both  explained  quantiutively  by  a  semiconductor 
theory  of  coherent  energy  transfer.  A  full  understanding  of  strong  injection  effects  on  the  semiconductor 
gain  curve  can  help  in  the  understanding  of  VCSEL  instabilities. 

The  VCSEL  lases  in  two  nearly  degenerate  cross-linearly  polarized  modes,  one  of  which  is 
stronger  by  a  factor  of  about  100  for  lasing  well  above  threshold.  A  weak  injection  may  switch  the 
polarization  of  the  VCSEL  lasing,  exchanging  the  roles  of  the  strong  and  weak  naodes.  Local  gain 
changes  induced  in  one  mode’s  spectrum  appear  in  the  VCSEL’s  other  polarized  mode  (but  not  in  higher- 
order  modes).  The  Rayleigh-gain  spatial  grating  written  in  one  polarization  is  also  seen  in  the  other 
polarization,  which  has  a  different  frequency  but  the  same  wavelength  and  wavevector  in  the  material; 
see  Fig.  3.* 

Subpicosecond  Nonlinear  Directional  Coupler* 

We  have  collaborated  with  Peyghambarian’s  group  to  demonstrate  subpicosecond  switching  and 
recovery  of  the  transmission  of  the  nonlinear  directional  coupler  with  cunent  injection,  which  we 
fobricated  previously  with  the  help  of  Boeing  Aerospace.  The  previous  performance  of  this  NLDC  with 
IB^s  pulses  was  very  little  better  than  that  of  a  passive  NLDC.  The  much  better  performance  with  lOQ- 
fi  pulses  occurs  at  the  no-gain/no-loss  point  (chemical  potential),  as  shown  earlier  in  a  single  waveguide 
laser  amplifier  by  Hultgren  and  Ippen.  The  physics  of  the  subpicosecond  recovery  is  still  controversial, 
and  we  and  others  continue  to  unravel  it. 


Figure  3.  The  no-injectioa  VCSEL  apectnun  (dished)  his  in  s-polirized  mode  and  i  p-polirized 
mode  which  ire  both  modified  by  the  p-poUrized  injection  (12  fiW  intncivity  power)  it  (thin 
line).  A  20  /iW  intncivity  injectkm  it  modifies  both  modes  and  produces  a  10  >tW  new  lasing 
line  It  Which  also  modifies  both  modes. 
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RESEARCH  FINDINGS 

Gaussian-Beam-Induced  Vortices  in  Vertical-Cavity  Surface-Emitting  Lasers*'^ 

We  have  triggered  transverse  pattern  modiflcatioDS  by  injecting  a  CW  monochromatic  (10  MHz) 
beam  into  a  vertical-cavity  surface-emitting  laser  (VCSEL).  The  injected  Gaussian  beam  has  a  frequency 
close  to  that  of  the  natural  TEMoo  mode  of  the  VCSEL  (Fig.  1).  The  spot  of  the  injection  is  not 
superimposed  with  the  original  lasing  spot  (Fig.  2).  However,  it  is  close  enough  to  lock  the  TEMn 
mode,  as  seen  in  the  interferogram  of  Fig.  3.  As  a  result,  we  observe  a  system  of  two  Gaussian  beams 
locked  in  frequency  and  in  phase,  propagating  on  parallel  axes,  and  with  waists  of  different  size  and 
location.  The  interferogram  of  Fig.  3  shows  two  vortices  of  opposite  topological  charges  located  on  both 
sides  of  the  system  of  two  Gaussian  beams. 

Rosanov  has  shown  theoretically  [Optika  i  Spektroskopiya  75,  861  (1993)]  that  if  such  a  system 
had  the  two  beams  propagating  along  the  same  axis,  a  series  of  vortices  would  appear,  with  their 
trajectories  along  concentric  rims  normal  to  the  axis  of  propagation.  These  vortices  do  not  display  the 
phase  singularity  behavior  in  die  plane  normal  to  the  axis  of  propagation  and  cannot  be  observed 
interferomecrically.  However,  when  die  beams’  axes  are  not  collinear,  vortices  appear  with  their 
trajectories  along  lines  oriented  in  the  same  direction  as  the  axes,  and  thus,  can  be  observed  as  phase 
singularities  in  the  plane  normal  to  die  axes  of  propagation.  In  our  experiment,  the  image  plane  coincided 
with  the  focal  plane  of  the  injected  beam,  at  die  VCSEL/air  plane.  In  this  plane,  in  a  bipolar  coordinate 
system  (r„  with  poles  coinciding  with  axes  of  the  beams,  the  locations  of  vortices  are  given  by  the 
following  formulas: 


r?  -  2  [l  ♦  [w  ♦  1]  a]  D*  ♦ 

L  I  DJ  [ 


,  m  -  0,  ±1,  ±2, 


rl  -  Wj 


_ ^ 

Wi  * 
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Figure  1 .  Injecting  (dashed  curve)  at  the  fieqtiency  of  the  TEMg,  mode  modifies  the  original  lasing  spectrum  (solid 
curve),  and  locks  the  laser  in  a  pure  TEMg,  mode  (dotted  curve). 

where  index  1  stands  for  the  laser  beam,  and  index  2  stands  for  the  injected  beam.  D  is  the  distance 
between  focal  planes  of  the  two  beams;  Wj  are  half-waists  of  the  beams;  Aj  are  beam  amplitudes;  and  X 
is  the  wavelength. 

The  pair  of  vortices  observed  in  the  experiment  corresponds  to  the  order  m  =  -IS.  A 
neighboring  pair  of  vortices  should  be  located  outward  at  a  distance  of  more  than  5  fim  from  the  observed 
vortices.  In  our  experimental  setup,  it  is  currently  out  of  observation.  The  distance  between  the  beams’ 
axes  is  not  a  parameter  in  the  above  formulas.  Hence,  bringing  the  beams  closer  increases  the  distance 
between  the  vortices,  and  vice  versa.  This  d^endence  was  clearly  conrirmed  in  the  experiment. 
Likewise,  rj  linearly  depends  on  w,.  When  w,  increases,  r,  increases  also,  separating  the  two  vortices 
further.  These  changes  in  positions  of  the  vortices  also  were  clearly  seen  in  the  experiment. 

We  conclude  that  our  VCSELs  do  not  emit  vortices.  The  vortices  observ^  with  injection  are 
an  interference  between  the  VCSEL  emission  and  the  injection.  Stationary  fringes  between  VCSEL  and 
injection-frequency  reference  beams  cannot  be  seen  without  injection  locking.  Furthermore,  efficient 
locking  requires  focusing  of  the  injection  onto  the  VCSEL.  For  these  reasons,  it  was  nontrivial  to 
conclude  that  the  vortices  arise  from  the  interference  of  two  beams  (VCSEL  and  injection)  with  waists 
in  the  same  plane,  rather  than  from  the  VCSEL  emission  itself.  This  is  another  example  of  the 
appearance  of  vortices  in  linear  optics,  just  as  in  optical  testing  experiments. 
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Figure  2.  The  original  lasing  spot  at  the  output  mirror  of  the  VCSEL  is  the  one  on  top  of  the  pattern.  The  injection 
has  been  positioned  underneath.  The  injection  is  close  enough  to  lock  the  VCSEL  lasing  in  fr^uency  and  in  phase. 


Figure  3.  An  interferogram  of  the  intensity  pattern  of  Fig.  2  shows  tw-o  vortices  of  opposite  topological  charge. 
These  vortices  move  up  and  outward  as  the  waist  of  the  injection  is  increased  or  as  the  injection  is  brought  closer 
to  the  original  lasing  spot  of  the  VCSEL. 
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SpRtiRl  SoUtons  in  r  Self*Foaising  Saniconductor  Gnin  Medium* 

A  beam  propagating  in  a  planar  waveguide  exhibits  the  characteristics  of  both  a  fundamental 
soliton  Oittle  change  in  waist  of  a  S-^m  beam;  diffraction  is  compensated  by  self-focusing  nonlinear 
refraction)  and  a  second-order  soliton  (10-/im  beam  focuses  to  3  /im  with  side  peaks).  These  one¬ 
dimensional  spatial  solitons  are  quite  stable;  bifurcations  of  transverse  solitary  waves  are  not  seen,  in 
sharp  contrast  with  the  case  for  two  transverse  dimensions.’*  Plasma-theory  planar-wave  beam- 
propagation  computations  agree  well  with  the  data  and  resemble  Kerr-medium  solitons  quite  closely. 


Compression  of  Bright  Opdcal  Pubes 

This  is  a  foeoretical  prediction  of  the  compression  of  bright  optical  pulses  by  copropagation  with 
a  dark  temporal  soliton,  an  extension  of  ideas  from  our  earlier  paper  on  i^uced  (spatial)  focusing  in  a 
semiconductor  defocusing  medium.’ 


Photon-Number  Correlations  Near  the  Threshold  of  Microcavity  Lasers*''* 

We  made  the  first  second-order  correlation  measurement  of  a  VCSEL  and  found  that  the  peak 
of  the  Fano-Mandei  parameter  K  ■  [(h*)  -  (n)’  -  (n)]/(n)  =  (n)[g'’'(0)  -  1]  occurs  at  the  same  pump 
intensity  that  the  derivative  of  the  output  intensity  becomes  discontinuous,  for  oat  0  ^  0.(X)dS  for  foe 
fraction  of  spontaneous  emission  going  into  foe  lasing  mode.  As  -•  1,  foe  discontinuity  in  output  power 
vanishes;  some  people  call  this  a  thresholdless  las^ .  We  suggest  foat  since  foe  low-pump  emission  is 
simply  spontaneous  emission,  it  is  better  described  as  LED  rafoer  than  laser  emission.  Using  two-level- 
atom  photon  statistics  wifo  GaAs  parameters,  we  show  foat  foe  intensity  fluctuations  of  foe  field, 
quantified  by  K,  still  exhibit  a  maximum  because  of  foe  transition  from  spontaneous  to  stimulated 
emission,  llie  transition  of  g^(0)  from  2  to  1  at  lasing  threshold  becomes  softer  and  softer  as  /?  -*  1. 
Therefore,  we  conclude  that  there  is  no  thresholdless  laser. 
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RESEARCH  FINDINGS 

In  this  project  we  proposed  to  perform  theoretical  and  experimental  studies  of  the  linear  and 
nonlinear  optical  properties  of  low-dimensional  semiconductor  structures,  such  as  quantum  dots,  quantum 
wells,  or  more  complicated  heterostructures.  Here  we  summarize  our  results. 

CdS  and  CdSe  Quantum  Dot  Research 

Nonlinear  absorption  spectra  of  quantum-confioed  CdSe  microcrystallites  in  glass  were  measured 
in  a  femtosecond  pump^robe  experiment  at  low  temperatures.  Simultaneous  bleaching  of  the  quantum- 
confined  transitions  and  an  induced  absorption  at  higher  energy  were  observed.  Spectral  hole  burning 
was  observed  as  the  pump  was  tuned  through  the  transition.  The  system  was  modeled  as  spherical 
crystallites  with  one  or  two  electron-hole  pairs  int^Kting  through  the  Coulomb  potential.  Our  numerical 
matrix  diagonal  ization  method  predicted  results  essentially  similar  to  the  experiment  and  attributed  them 
to  state  filling  of  the  one-pair  states  and  generation  of  two-pair  states. 

The  linear  absorption  spectra  of  three  samples  of  CdS  quantum  dots  in  glass  at  10  K  are  shown 
in  Fig.  1.  The  spectrum  labeled  "bulk"  refers  to  a  glass  wiA  semiconductor  microcrystallites  large 
enough  to  retain  the  three-dimensional  bulk  properties.  This  spectrum  is  typical  of  a  bulk  CdS  absorption 
spectrum  with  a  sharp  exciton  and  structureless  Coulomb-enhanced  continuum  absorption  at  higher  photon 
energies.  The  average  crystallite  size  becomes  smaller  for  lower  heat-treatment  temperatures.  The 
spectra  labeled  'QD*  refers  to  samples  with  small  crystal  sizes.  The  quantum-confinement  effects  are 
clearly  observable  in  this  sample.  The  absorption  has  shifted  to  higher  energies,  as  expeaed  for  the 
confinement  effect.  Furthermore,  discrete  quantum-confined  electron-hole-pair  states  appear  in  the  "QD* 
sample. 

An  example  of  measured  optical  nonlinearities  for  quantum  dots  in  glass  is  shown  in  Fig.  2.  The 
absorption  change  for  resonant  excitation  into  die  energetically  lowest  quantum-confined  transition  was 
measured.  The  top  spectrum  in  Fig.  2  displays  the  linear  absorption  of  a  CdS-doped  glass  quantum-dot 
sample  at  T  =  10  K.  This  spectrum  has  two  peaks  at  energies  of  approximately  2.9S  eV  a^  3.26  eV. 
The  low-energy  peak  originates  from  the  transition  between  the  Is-hole  and  Is-electron  states.  The  high- 
energy  peak  is  due  to  the  transition  between  the  Ip-hole  and  Ip-electron  states. 


Figuiel.  The  linear  absorptioo  spectn  at  T  a  10  K  of  CdS  crystallites  in  glass.  The  two  spectra 
labeled  "QD*  refer  to  two  samples  with  small  crystal  sizes  in  the  qusntum-confinenieot  regime. 
The  *bulk”  sample  has  large  crystal  sizes,  exhibiting  bulk  CdS  properties.  (The  samples  were 
prepared  by  Liu  and  Risbud.) 
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Figure  2.  The  top  curve  is  the  Uneer  sbeorptiou  qMctium  of  a  sample  of  CdS  quantum  dots  in 
glass.  The  bottom  trace  shows  the  energetic  position  of  the  pump  pulse.  The  intermediate  trace 
is  the  change  in  the  linear  absorption  spectrum  as  a  result  of  excitation  by  the  pump  pulse,  d  is 
the  thickness  of  the  sample.  (Tte  sample  was  prepared  by  Liu  and  Ridnid.) 


The  narrow  peak  at  the  bottom  of  Fig.  2  displays  the  spectral  position  of  die  pump  pulse.  A 
broadband  probe  pulse  detects  the  changes  made  in  the  absorption  spectra  as  a  result  of  excitation  by  die 
pump  pulse.  Probe  absorption  without  previous  excitation  by  the  pump  is  equivalent  to  the  generation 
of  one  electron*hole  pair  in  the  unexcited  quantum  dot.  If,  however,  a  pump  photon  has  been  absorbed 
previously,  the  probe  beam  generates  a  second  electrDn*faole  pair  in  the  presence  of  the  pump-generated 
pair. 

The  lower  spectrum  in  Fig.  2  shows  die  measured  changes  in  die  absorption,  -Aod  =  d[o^ 
(without  pump)  -  (with  pump)].  In  such  a  spectrum,  a  positive  peak  corresponds  to  bleaching  and 
a  negative  p^  i^icates  an  induced  absorption,  i.e.,  an  increase  in  the  probe  absorption  as  a 
consequence  of  die  presence  of  the  pump-generated  electron-hole  pair.  The  nonlinear  spectrum  in  Fig.  2 
shows  a  positive  pe^  around  the  Is^  -  Is.  transition,  centered  around  die  pump,  and  a  n^ative  peak  on 
the  high-energy  side,  centered  around  3.178  eV.  The  bleaching  of  the  Is^  -  Is,  transition  is  the  result 
of  state  fillmg.  The  generation  of  one  electron-hole  pair  by  the  pump  causes  saturation  of  the  one-pair 
transition  (bleaching  of  the  linear  absorption  peak). 

The  origin  of  the  induced  absorption  feature  (the  negative  peak  in  die  lower  spectrum  of  Fig.  2) 
is  assigned  to  the  generation  of  two^iair  states  in  the  quantum  dot.  The  pairs  are  creked  by  absorption 
of  one  pump  and  one  probe  photon.  The  theoretical  analysis  of  the  linear  and  third-order  nonlinear 
optical  properties  has  to  include  the  one-  and  two-electron-bole-pair  states.  These  investigations 
consistendy  lead  to  the  conclusion  that  Coulomb  effects  are  important  even  for  the  smallest  quantum  dots. 
Using  numerical  matrix  diagonal  ization  techniques,  the  energies  and  wave  functions  were  obtained  for 
the  one-  and  two-pair  ground  states  and  for  all  the  excited  pair  states.  With  these  wave  functions,  the 
various  dipole  matrix  elemmits  were  evaluated  for  transitions  between  die  ground  state  and  the  one-  and 
two-electron-hole-pair  states.  The  changes  in  the  absorption,  -Aod,  were  then  calculated,  as  shown  in 
Fig.  3.  Like  the  experiment,  the  theory  also  exhibits  a  decreasing  absorption  feature  (bleaching)  around 
the  lowest  quantum-confined  transition  and  an  increasing  absorption  feature  on  die  high-energy  side. 
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Figure  3.  Computed  ebeotptioa  changes  for  a  semiconductor  quantum  dot  of  radius  R  »  a^  (where 
a,  is  the  exciton  Bohr  radius),  assuming  a  pump  and  probe  geometry,  and  pumping  into  the 
energetically  lowest  one-pair  state. 

Optical  Nonlincarities  of  CdTe  Quintum 

We  extended  our  quantum  dot  research  to  include  microcrystallites  of  CdTe  in  glass  matrices. 
CdTe  has  the  largest  Bohr  radius  of  any  semiconductor  successfully  grown  as  quantum  dots  in  glass. 
Therefore,  quantum-confinement  effects  can  even  be  seen  in  larger  ^ts. 

To  distinguish  between  intrinsic  effects  of  photogenerated  electron-hole  pairs  and  effects  of 
trapped  charges,  we  measured  the  nonlinear  absorption  changes  at  early  times  using  femtosecond  laser 
pulses.  Femtosecond  pulses  at  620  nm  were  generated  in  a  colliding^pulse  mode-locked  dye  laser 
cavity  and  were  amplified  by  an  8.5  kHz  copp«  vapor  laser.  Figure  4(a)  shows  the  linear  abwrption 
spectrum  of  one  of  the  CdTe  quantum  doc  samples  with  average  crystallite  size  of  36  A.  The  Sa  a^  Sg 
transitions  in  diis  figure  correspond  to  two  quantum-confined  transitions  that  share  the  same  electron 
levels.  The  measured  absorption  (5aL)  changes  immediately  aftw  excitation  inside  the  I'me,  and  the 
pump  spectral  position  is  displayed  in  Fig.  4(b).  The  boL  spectrum  consists  of  two  features 
corresponding  to  bleaching  of  both  Sa  and  lines.  The  bleached  Sb  line  has  larger  magnitude  compared 
with  the  Sa  bleachmg,  in  spite  of  the  fact  that  the  linear  Sa  absorption  is  larger  than  Sb  [see  Fig.  4(a)]. 
These  observations  are  reproduced  by  a  sinqile  model  whose  result  is  presented  as  the  da^ed  line  in  foe 
inset  of  Fig.  4.  In  this  modd,  we  take  into  account  transitions  between  foe  ground  state  (no  dectron-hole 
pair)  and  foe  Sa  and  Sb  one-pair  states.  In  addition,  we  also  include  transitioos  to  foe  excited  two-pair 
states,  which  are  energetically  above  foe  energy  of  twice  foe  one-pair  transitions.  The  size  distribution 
of  quantum  dots  is  modded  as  an  inhomogeneous  broadening  of  all  foe  transitions.  The  qualitative 
agreement  between  Fig.  4(b)  and  the  dashed  spectrum  of  foe-  inset  shows  that  state  filling  of  one-pair 
transitions  is  clearly  responsible  for  foe  observed  absorption  saturation. 

To  further  examine  foe  state-filling  effect,  we  tuned  foe  pump  pulse  inside  foe  Sa  transition,  as 
shown  in  Fig.  4(c).  As  expected,  we  see  an  Sa  to  Sb  bleaching  ratio  that  is  not  only  larger  than  that  of 
Fig.  4(b),  but  also  larger  tl^  foe  linear  Sa  to  Sb  absorption  coefficient  ratio.  Furthermore,  an  induced 
absorption  feature  is  apparent  on  foe  high-energy  side  of  foe  Sb  bleaching.  As  shown  by  our  model 
calculation  in  foe  inset  of  Fig.  4  (foe  solid  curve),  this  induced  absorption  is  foe  result  of  foe  two- 
electron-hole-pair  transitions.  The  larger  bleaching  of  Sa  compared  with  Sb,  and  foe  induced  absorption 
on  foe  high-energy  side  in  foe  spectrum  shown  by  the  solid  line  of  foe  inset,  are  clearly  similar  to  foe 
experiment  of  Fig.  4(c). 
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Figure  4.  («)  T  inrur  abaoiptioa  of  the  CdTe  <]uaBtum-dot  aaaiple  at  IS  K.  (b)  and  (c)  Meaaured 
abaorptiao  chaagea  for  puiD|HOg  at  at  620  and  »  653  fun,  reapactivdy.  Tba  apactnl  focatkm  of 
the  pompe  are  ahown  oa  the  upper  pert  of  each  figure.  Inaet:  cafculated  afaaorptioa  changea  for 
the  two  pumping  aituatioaa.  The  pump  poaitioaa  are  diaplayed  by  arrows.  The  dashed  spectrum 
(solid  spectrum)  is  calculated  for  a  detming  of  6  (zero),  suitable  for  compariaoa  with  (b)  ((c)]. 

(The  aam|de  was  prepared  by  Liu  and  Riabud.) 

We  proposed  the  following  scenario  to  explain  the  observed  phenomena.  The  first  photon  creates 
an  electron-^ole  pair  inside  the  quantum  dot.  The  pair  creation  may  happen  with  or  widxMit  die  presence 
of  charged  impurities  or  trap  states.  Since  the  dectron  and  hole  interact  via  the  attractive  Coulomb 
potential,  their  wave  fonctioos  deviate  significantly  from  those  of  the  strong-confinement  ^roximation 
(no  Coulomb  interaction).  There  is  strong  evidence  that  the  long-time  behavior  of  die  absorption  changes 
are  determined  by  carriers  trapped  at  surfKe  states  (e.g.,  Cd’^  ions)  or  into  localized  states  in  die  glass 
outside  die  quantum  dot.  The  corre^nding  absorption  changes  result  from  die  Cfoulomb  imeraction  of 
the  probe-goierated  electron-hole  pairs  with  the  trapped  charges  and  state  filling,  whidi  may  become 
partial  if  a  carrier  leaves  die  dot.  Therefore,  state  filling  and  Coulomb  interaction  are  responsible  for 
absorption  dianges  at  early  times  after  excitation,  and  Coulomb  effects  are  the  only  important  mechanism 
for  later  times.  Both  photoinduced  absorption  changes  and  photoluminescence  of  traps  have  long-lasting 
conqiooents  yidding  evidence  diat  long-lived  nonlinear  optical  conqionents  are  associated  with  traps. 


Femtosecond  Dynamics  of  Type  II  Quantum  Wdls 

The  many-body  effects  associated  widi  a  two-component  plasma,  i.e.,  dectrons  and  holes,  have 
been  studied.  In  type-n  quantum  wdls,  on  the  other  hand,  the  dectrons  are  removed  from  the  holes,  and 
the  dynamics  of  holes  or  dectrons  can  be  investigated  separatdy.  For  exanqile,  whm  the  GaAs  layer 
tiiickness  is  less  dian  35  A  and  AlAs  layer  thickness  is  greater  than  16  A,  the  lowest  T-dectron  state  of 
tiie  GaAs  wdl  lies  energetically  above  the  X-minimum  of  the  AlAs  barrier,  whereas  the  lowest  energy 
hole  state  is  still  in  die  GaAs  layer  (see  the  inset  of  Fig.  5).  As  a  result  of  such  a  band  dignment,  optical 
excitation  of  dectron-hole  pairs  in  the  GaAs  layer  is  followed  by  a  spatid  separation  of  die  two  plasma 
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components;  electrons  scatter  from  the  T-state  in  GaAs  to  the  X-sute  in  AiAs.  This  allows  the 
measurement  of  the  dynamics  of  the  hole  relaxation  in  the  GaAs  layer  where  they  are  spatially  separated 
from  the  electrons  in  the  AlAs  layer.  We  have  made  such  measurements  at  T  »  10  K  using  130-f$  laser 
pulses,  tuned  resonantly  inside  the  heavy-hole  exciton  peak  of  a  type-U  quantum-well  structure  consisting 
of  a>28-A-thick  GaAs  (10  monolayers)  and  a57-A-thick  AlAs  (20  monolayers)  layers. 

Figure  S  displays  the  absorption  spectra  of  the  sample  at  various  time  delays  between  the  pump 
and  probe  pulses  together  with  the  pump  spectrum.  A  bleaching  and  blue  shift  of  the  absorption  spectrum 
at  both  the  light-hole  and  the  heavy-hole  excitons  is  clearly  observed  for  the  zero-ps  time  delay.  Blocking 
of  the  conduction  and  heavy-hole  states  by  the  electrons,  which  have  not  yet  left  the  GaAs  layer,  is 
responsible  for  this  effect.  In  the  3-ps  spectrum,  die  li^t-hole  exciton  blue  shift  is  completely  recovered, 
while  die  heavy  hole  is  still  bleached  and  shifted.  This  behavior  is  caused  by  the  faa  that  the  F-X 
electron  scattering  has  already  taken  place  and  only  the  F-point  heavy-hole  phase-space  filling  is  left. 
The  relaxation  (cooling)  of  the  heavy  holes  contributes  to  Ae  behavior  of  the  l(X)-ps  trace.  Also,  we 
observe  a  complete  recovery  of  the  light-hole  exciton  peak  due  to  decreased  screening  of  'cool'  heavy 
holes.  As  the  holes  relax  to  the  lowest  energy  states,  the  bleaching  of  the  spectral  region  on  the  low- 
energy  side  of  the  heavy-hole  exciton  (see  the  region  around  X  =>  670  nm)  is  increased,  while  the 
bleaching  of  the  high-energy  side  of  the  heavy-hole  exciton  (see  the  region  around  X  ^  660  nm)  is 
recovered. 

In  our  theoretical  analysis  of  the  many-body  effects  in  the  band-edge  absorption  spectra  of  highly 
excited  type-I  and  type-U  semiconductor  quantum-well  structures,  we  assumed  perfect  electron-hole 
charge  separation,  so  that  only  the  F-point  holes  contribute  to  phase-space  filling  of  the  GaAs  exciton 
states.  We  numerically  solved  the  mterband  polarization  equation,  including  inhomogeneous  broadening 
due  to  well-width  fluctuations  by  averaging  ^e  spectra  over  a  distribution  of  well  thicknesses. 


Figure  S.  Absoiptioa  spectra  at  differeot  time  delaya  for  reaooaiit  heavy-hole  exciton  pumping 
in  a  GaAa-AlAa  type-U  quantum  well.  The  inaei  ihoara  die  energy  level  diagram  for  the  sample. 
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Figure  6  shows  computed  quasi-equilibrium  absorption  spectra  for  different  plasma  densities.  In 
the  low-temperature  spectra  of  Fig.  6(a).  saturation  and  blue  shifting  of  the  HH  exciton  and  a  slight  red 
shift  of  the  LH  exciton  can  be  seen.  For  elevated  plasma  temperatures.  Fig.  6(b)  shows  that  the  exciton 
blue  shift  is  substantially  reduced,  since  the  holes’  distributions  involve  more  band  states,  substantially 
reducing  Coulomb  enhancement  and  Pauli  blocking  effects.  The  type-I  results  in  Fig.  6(c)  also  show  HH 
exciton  saturation  and  blue  shift,  along  with  a  blue  shift  of  the  LH  exciton  due  to  phase-space  filling  by 
die  electrons. 


Figure  6.  Computed  absorpdoo  qiectn  for  type-I  ud  type-II  quantum  wells  with  GaAs  well 
thickness  of  30  A.  The  spectra  are  for  the  plasma  densities  na^o  «  0.2  (short  dashed),  0.4 
(inedium  dashed),  and  0.8  Ooog-short  dashed).  The  linear  spectra  are  shown  as  solid  lines,  and 
So  is  the  bulk-exciton  Bohr  radius,  (a)  Type-Il,  plasma  temperature  T  »  30  K;  (b)  Type-II,  T  = 
100  K;  (c)  Type-I,  T  -  30  K. 
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There  is  do  gain  in  the  type-II  spectra,  whereas  clear  regions  of  negative  absorption  develop  in 
the  type*I  spectra.  This  absence  of  gain  in  moderately  excited  ideal  type-II  structures  is  a  direct 
consequence  of  the  spatial-electron  and  hole-plasma  separation.  The  gain  part  of  a  type-I  spectrum  is 
replaced  by  a  zero-absorption  region  in  type-II  structures  where  the  quasi-chemical  potential  coincides 
with  the  onset  of  absorption  for  sufficiently  high  densities. 

Photon  Echo  in  Semiconductors 

We  began  the  theoretical  study  of  the  photon  echo  in  semiconductors  using  the  semiconductor 
Bloch  equations.  For  the  photon  echo  case,  we  assumed  a  configuration  where  the  two  exciting  pulses 
are  incident  under  an  angle.  We  took  the  first  pulse  to  be  weak  so  that  we  could  assume  linear  response. 
The  strong  second  pulse  was  treated  in  all  orders.  The  center  frequencies  of  both  pulses  were  at  the  Is 
exciton  resonance,  assuming  pulse  separation  ofr  -  400  fs.  Figure  7  shows  that  for  the  case  of  a  weak 
first  pulse,  an  almost  instantaneous  signal  and  no  photon  echo  at  -i-400  ps  occurs.  Our  analytical 
calculations  proved  that  this  signal  is  due  solely  to  the  exchange  correlation  b^een  the  excited  excitons. 
For  higher  pulse  intensities,  we  see  the  gradual  development  of  an  echo  signal,  which  coexists  with  the 
instantaneous  signal  for  intermediate  intensities. 

To  analyze  the  origin  of  this  scenario.  Fig.  8  plots  the  time  dependence  of  the  renormalized 
bandgap  for  the  excitation  conditions  of  Fig.  7.  Compar'mg  Figs.  7  and  8  reveals  that  the  echo 
contribution  in  the  time-resolved  signal  occurs  as  soon  as  the  continuum  states  are  shifted  into  resonance 
during  the  presence  of  the  first  pulse  (bandgap  shift  below  -lEn  in  Fig.  8).  Consequently,  direct 
continuum  excitation  is  possible,  which  yields  a  photon  echo  signal  at  4(X)  fs  because  of  the  intrinsic 
inhomogeneous  broadening  of  the  electron-hole  continuum  sutes. 
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Figure  7.  Time-reeolved  signet  in  the  photon-echo  direction  for  increasing  strength  of  the  first 
pulse,  E),  Excitation  occurs  at  the  exciton  resonance.  The  time  200  fs,  the  time  delay 

r  400  fi,  and  the  pulse  FWHM  is  100  fi  for  both  pulses.  The  peak  value  of  the  dipole 
coupling  energy  of  the  second  pulse  is  d^  E|  ~  0. 1  E^,  where  the  exciton  binding  energy  Et  > 
16  meV  in  Cd^.  The  correqnnding  peak  amplitudes  of  the  first  pulse  are  shown  with  increasing 
dash  length  for  dw  E,  »  0.01  E^,  0.03  E^,  0.05  E^,  0.07  Et.  and  for  0.1  Et  (solid  line). 


Figure  8.  Renormalized  band  edge  as  a  ftmctioo  time.  E*o  is  the  unrenormalized  band  edge.  The 
different  curvea  are  d^  E,  »  0.01  Et,  0.03  Et,  O.OS  Et,  0.07  Et,  and  0.01  Et,  from  top  to 
bottom,  respectively.  All  parameters  are  the  same  as  in  Fig.  7. 
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RESEARCH  FINDINGS 

The  results  we  obtained  during  the  course  of  this  study  are  summarized  below. 

Adiabatic  Following  in  Semiconductors 

Adiabatic  following  is  an  off-resonant  effect  which  occurs  when  the  duration  of  a  light  pulse  is 
less  than  the  phase  relaxation  time  T,  and  the  magnitude  of  the  pulse  detuning  is  greater  than  its 
inhomogeneous  linewidth.  Under  diese  conditions,  the  quantities  in  the  optical  Bloch  equations  which 
describe  the  system,  namely  the  inversion  and  polarizatioo,  have  a  time  dependence  detmnined  by  the 
instantaneous  amplitude  of  the  light^iulse  envelope,  i  .e..  they  follow  the  fidd. 

The  above  conditions  on  the  pulse  duration  and  detuning  can  be  satisfied  in  a  semiconductor  using 
femtosecond  laser  systems.  This  is  most  easily  achieved  for  the  exciton  resonance  in  GaAs  multiple- 
quantum-well  structures  (MQWs),  which  has  a  coherence  time  approaching  one  picosecond.  The 
response  of  the  exciton  to  a  nonresonant  light  pulse  can  be  describ^  by  the  inversion  and  ihterband 
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polarization  of  die  system  using  the  semiconductor  Bloch  equations,  which  include  the  many-body 
Coulomb  effects  in  time^ependent  Hartree-Fock  approximation.  The  numerical  solution  of  these  coupled 
equations  not  only  explains  the  Stark  shift,  as  shown  previously,  but  also  indicates  that  the  system’s 
inversion  adiabatically  follows  the  field.  This  aspect  of  the  Stark  effect  is  manifested  experimentally  in 
the  time-resolved  absorption  measurements  as  a  fast  bleaching  recovery  of  the  exciton  line. 

Time-resolved  pump-probe  experiments  were  carried  out  using  a  synchronously  pumped  mode- 
locked  dye  laser  with  an  average  output  power  of  2S  mW,  82  MHz  repetition  rate,  and  center  wavelength 
tunable  ftom  850  to  870  nm  for  room-temperature  experiments.  A  colliding  pulse  mode-locked  (CPM) 
dye  laser  amplified  by  copper  vapor  lasers  (CVL)  in  cascade  operating  in  750  nm  to  800  nm  was 
employed  for  low-temperature  (a>  10  K)  experiments.  For  this  purpose,  we  first  generated  a  continuum 
using  CPM  pulses  amplified  by  CVL,  and  then  the  near-IR  pulses  were  obtained  by  reamplifying  the 
desired  portion  of  the  generated  continuum  using  a  second  CVL.  The  autocorrelation  of  the  pump  pulse 
and  the  cross  correlation  of  the  pump  and  probe  pulses  in  this  case  were  200  fs  and  300  ft,  respectively. 
The  samples  were  molecular  beam  epitaxy  (MBE)  grown  GaAs/Al,Ga,.^s  multiple-quantum-well  (MQW) 
and  multiple-coupled-quantum-well  (MCQW)  structures  with  various  well  widths  and  barriers.  In  all 
room-temperature  measurements  the  laser  center  wavelength  was  adjusted  to  allow  detuning  between  four 
and  five  ^  below  the  heavy-hole  exciton.  Here  Er  is  the  bulk  GaAs  Rydberg  energy,  i.e.,  4.2  meV. 
The  detunings  for  low-temperature  measurements  were  smaller.  The  spectral  transmission  of  the  probe 
through  the  sample  for  different  time  delays  between  pump  and  probe  was  measured  by  an  optical 
multichannel  analyzer  at  the  output  of  a  spectrometer. 

Figure  1(a)  shows  the  low-temperature  absorption  spectrum  of  a  GaAs  MCQW  sample  for 
pumping  below  the  exciton  resonance  a^  for  different  time  delays,  tp,  where  t,  =  t(probe)  -  t(pump). 
The  solid  curve  represents  the  linear  absorption,  while  the  dotted  spectrum  corresponds  to  t,  0.  It 
clearly  shows  that  the  heavy-hole  exciton  has  been  blue  shifted  and  bleached.  The  dashed  curve  in 
Fig.  1(a),  which  corresponds  to  ^  =  500  ft,  demonstrates  that  the  blue  shift  and  bleaching  are  mostly 
recovered.  The  complete  recovery  takes  nanoseconds  as  a  result  of  carrier  generation  caused  mainly  by 
the  spectral  overlap  of  the  pump  the  san^ile’s  absorption  spectrum.  The  transient  exciton  blue  shift 
is  the  manifestation  of  the  opti<^  Stark  effect.  This  bleaching  recovery  is  the  signature  of  the  transient 
adiabatic  following. 

For  the  analysis  of  our  experiments,  we  used  the  semiconductor  Bloch  equations,  i.e.,  the  coupled 
equations  of  motion  of  the  expectation  value  of  the  population  of  the  state  k  and  the  interband 
polarization. 

We  solved  foe  semiconductor  Bloch  equations  numerically  for  different  time  delays,  assuming 
constant  dephasing  and  carrier-relaxation  rates.  The  calculated  spectrum  for  comparison  with  the  room- 
temperature  data  of  Fig.  1(a)  is  shown  in  Fig.  1(b).  As  in  the  experiment,  foe  pump-pulse  duration  and 
detuning  satisfied  foe  adiabatic  following  conditions  mentioned  earlier.  The  temporal  behavior  of  the 
exciton  in  Fig.  1(b)  shows  good  qualitative  agreement  with  foe  data.  That  is,  the  exciton  both  bleaches 
and  shifts  at  negative  time  delays.  The  Stark  shift,  which  reaches  a  maximum  at  a  negative  time,  fully 
recovers  after  several  hundred  femtoseconds,  while  the  bleaching,  which  is  maximized  at  t,  =  0,  does 
not  quite  completely  recover. 

Figure  2  shows  foe  temporal  behavior  of  foe  created  carrier  density  along  with  the  pump-pulse 
intensity.  Note  that  similar  to  foe  experiments,  there  is  a  fast  component  of  foe  density,  which  follows 
foe  pump,  and  a  small  long-lasting  tail  due  to  the  incoherent  component  of  foe  real  carrier  generation. 
The  fast  conqwnent  is  a  result  of  foe  coherent  response  of  the  carrier  density.  It  is  foe  transient  presence 
of  this  density  which  is  responsible  for  foe  fast  bleaching  recovery.  We  assign  this  behavior  to  the 
ultrafast  adiabatic  following  in  semiconductors. 


37 


CAUUER  KELAXATION  STUDIES... 

pbyghambarian 


-7  -5  -3  -1  1 

Ensrgy  (E-E^/Ei 

(b) 


Fiiura  1.  (■)  McMorad  abntptioB  apectn  of  a  GaAs-AlGaAs  MCQW  at  T  -  10  K  at  difimt 
tUM  dakjn.  (b)  Caloilated  abaoiptioo  qwctta  for  coiapariaoa  with  Fig.  1(a). 
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Figure  2.  Calculated  temporal  behavior  of  the  carrier  density  generated  by  pun^  pulse. 


Femtosecond  Nonequilibrium  Ctnier  Reltxatioa  in  Bulk  CdSe 

A  femtosecond  spectral  hole-burning  technique  was  employed  to  study  the  relaxation  of 
nonequilibrium  carriers  via  carrier-carrier  (CCS)  and  carrier-LO  phonon  (CPS)  scattering.  Excitation 
by  70-6  laser  pulses  several  LO  phonon  energies  above  the  exciton  resonance  at  10  K  in  CdSe  results 
in  a  transient  spectral  hole  that  dis^pears  in  less  than  100  fs.  At  the  onset  of  the  pump  pulse  we 
observed  a  nondiermal  distribution  dik  essentially  extended  from  the  pump  energy  of  1.99  eV  to  lower 
energies,  indicating  participation  of  both  CCS  and  CPS  as  expected.  The  eiqieriments  were  performed 
with  70-£i  pulses  generated  at  1.99  eV  from  an  amplified  colliding  pulse  mode-locked  dye  laser. 
Figure  3(a)  shows  the  linear  absorption  spectrum  of  the  sample  at  10  K.  The  two  excitonic  peaks  labeled 
A  and  B  originate  from  the  heavy-bole  and  light-hole  valence  bands  split  by  crystal  field  interaction. 
Absorption  changes,  -Aa,  following  excitation  by  the  pump  pulse  were  measured  as  a  function  of  time 
delay  between  the  punq>  and  probe  pulses.  Figure  4(a)  shows  -Aa  spectra  in  50-fs  intervals.  The  0-fs 
and  SB-fs  spectra  show  the  presence  of  a  hot  nonfoermal  spectral  hole,  as  indicated  by  die  batdied  area, 
on  the  high-energy  side  of  A  and  B  exciton  bleaching.  As  die  nonthermal  distribution  thermalizes,  the 
spectral  hole  washes  out  and  only  die  bleached  excitons  remain. 

Our  dieoretical  analysis  was  based  on  the  semiconductor  Bloch  equations.  CCS  and  CPS  were 
included  in  die  carrier  collision  rates,  and  screening  was  treated  quasi-statically.  Results  of  our  theory 
are  presented  in  Figs.  3(b)  and  4(b).  Figure  3(b)  shows  the  cidcuiated  linear  absorption  spectrum, 
consisting  of  die  A  and  B  excitons.  The  0-fi  spectrum  in  Fig.  4(b)  clearly  displays  the  spectral  hole 
peaked  at  die  punqi  position  12  E^,  where  Ek  =  1S.7S  meV,  the  A-exciton  Rydberg  energy.  The  two 
excitons  are  conqilie^y  bleached.  The  hole  has  a  tail  extending  to  low  energies  with  respea  to  pump 
energy,  as  was  (^served  in  Fig.  4(a),  in  good  agreement  with  the  experiment.  The  theory  also  ^ows 
that  the  spectral  hole  washes  out  by  100  fi,  and  the  excitons  suy  bleached  for  later  times,  indicating  very 
large  carrier-carrier  scattering  rates.  Our  analysis  shows  how  the  very  short  scattering  times  (both  the 
CCS  and  the  CPS)  yidd  an  almost  instantaneous  bleaching  of  both  die  exciton  resonance  and  the 
Coulomb  enhancement  of  die  lower  continuum  states.  Since  the  LOiihonon  emission  processes  are 
essentially  as  fist  as  the  CCS,  the  Pauli  blocking  (spectral  hole)  primarily  affects  the  low-energy  side  of 
the  pun^)  frequency.  After  iqiproximatdy  a  hundr^  femtoseconds,  the  phase-space  blocking  is  similar 
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Figure  3.  (.a)  The  meewied  lineer  abaoiptioa  of  our  CdSe  platelet  at  10  K.  (b)  The  calculated 
linear  aboerpdon. 


PHOTON  ENERGY  («V)  ENERGY  (Ru  -  E,)/E« 

Figure  4.  (a)  The  measured  diaage  in  the  abeotpoon  coefRcient  observed  under  the  1.99  eV 
excitatkm  at  10  K.  The  time  delay  of  each  spectnim  is  ihown  in  the  figure.  The  hatched  area 
shows  the  nonthermalized  distribution  of  hot  earner*,  (b)  Calculated  pump*induced  absorption 
for  various  time  delays  bewteen  die  pusqi  and  probe  pwlaaa. 
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to  that  of  a  high-temperature  thermal  plasma  where,  in  this  case,  the  temperature  range  is  essentially 
given  by  the  LO-phonon  energy.  This  is  a  result  of  both  CPS  and  CCS,  where  the  CPS  reduces  the  mean 
kinetic  energy  of  the  plasma,  and  the  CCS  yields  a  quasi-thermal  distribution  where  even  the  states  below 
Awuo  ve  filled.  The  most  prominent  effect  of  the  phase-space  blocking  of  the  low  k-sutes  is  the 
reduction  of  the  absorption  increase  below  the  A-«xciton,  which  in  the  experimental  data  sets  in  after 
100  fs.  In  this  investigation,  the  experimental  results  seem  to  be  well-explained  within  the  simple  model 
of  equilibrium  phonon  distributions.  An  enhanced  temporal  resolution,  however,  might  prove  the 
signiftcance  of  oonequilibrium  and/or  coherent  phonon  effects. 


Coherent  Phonon  ExdUtion  in  Semiconductors 

Bil,  is  a  layered  semiconductor  with  two-dimensional  excitons  trapped  by  sucking  faults.  These 
excitons,  labeled  as  R,  S,  and  T,  are  extremely  sharp,  with  long  Tj  times,  as  shown  in  Fig.  S. 
Monitoring  the  absorption  change  at  the  peaks  of  R,  S,  and  T  excitons  resulted  in  an  interesting 
observation  shown  in  Fig.  6.  An  oscillatory  signal  is  detected  at  the  R,  S,  and  T  peaks  and  at  a 
wavelength  of 630  nm,  which  is  in  the  transparency  region  of  the  semiconductor.  We  verified  that  these 
oscillations  are  not  the  result  of  Rabi  flopping  bec^e  of  the  presence  of  the  signal  at  630  nm,  and  also 
because  the  signal  lasts  for  ac30  ps.  They  cannot  be  the  result  of  quantum  beats  between  the  excitons 
either  because  the  frequency  of  oscillations  does  not  correspond  to  the  energy  s^aration  of  any  of  the 
excitons.  The  energy  separation  of  a  14.3  meV  observed  in  these  oscUIations  precisely  matches  the  TO 
phonon  frequency  in  Bil,,  Thus,  we  concluded  that  our  dau  arise  from  a  coherent  excitation  of  TO 
phonons  in  the  medium,  presumably  by  impulse  stimulated  Raman  scattering  process.  The  »  30-ps  decay 
time  for  the  oscillations  then  is  due  to  the  vibrational  dephasing  as  a  result  of  coupling  with  incoherent 
phonons  in  this  two-dimensional  system. 


Figures.  Linear  abeorptioaspectnun  of  Bil,  at  T  >  10  K.  R,  S,  and  T  correapood  to  the  three- 
stacking-fault  excitons  in  this  layered  semiconductor. 
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Figure  6.  The  meesured  change  in  abaoq)don  for  the  R,  S,  and  T  excitona  and  the  transparency 
region  at  630  nm  of  the  Bill  sample  of  Fig.  5. 

Femtosecond  Pulse  Propagation  in  Semiconductor  Waveguides 

We  observed  a  new,  interesting  effect  that  we  called  'propagation-induced  escape  from  adiabatic 
following”  in  a  GaAs-AlGaAs  waveguide.  This  observation,  which  agreed  wdl  with  our  theory,  indicates 
that  high-intensity  femtosecond  pulses  break  up  as  they  propagate  through  a  semiconductor  when  they  are 
tuned  below  the  exciton  resonance.  This  pdse  breakup  results  from  the  escape  of  the  system  from 
adiabatic  following  regime  as  the  pulse  propagates  throu^  the  material.  The  summary  of  results  of  this 
work  may  be  found  in  P.  Harten  et  al.,  Phys.  Rev.  Lett.  69,  852  (1992). 


Spectral  Hole  Burning  in  the  Gain  of  an  Inverted  Semiconductor 

We  have  obtained  unambiguous  experimental  evidence  for  the  presence  of  a  spectral  hole  in  the 
gain  region  of  an  inverted  semiconductor.  A  three-beam  femtosecond  experiment  was  performed  for  this 
study.  The  first  femtosecond  laser  beam  was  used  as  the  ’gain”  beam,  whose  function  was  to  gen^ate 
high-density  electron-hole  pairs  to  invert  the  semiconductor  and  produce  the  gain  region.  Thus,  the 
material  was  pumped  optically,  not  electrically.  The  two  other  femtosecond  laser  beams  were  the  pump 
pulse  to  punq>  inside  the  gain  r^ion  and  a  broadband  probe  pulse  to  detect  foe  changes  made  by  foe 
pump.  The  broad  probe  allowed  monitoring  of  foe  entire  gain  region  of  foe  semiconductor 
simultaneously.  We  cfoserved  spectral  boles  that  persisted  for  foe  duration  of  foe  pump.  These  holes 
moved  with  foe  pump  frequency  as  foe  pump  tuned  across  foe  gain  r^ion,  as  expected.  Figure  7  shows 
a  typical  example  of  foe  hole-burning  results  where  foe  presence  of  foe  burned  hole  is  clearly  seen,  and 
its  distq>pearance  in  less  foan  200  fs  is  denoonstrated.  llie  incident  pulse  duration  was  » 100  fs. 


42 


CAUUEft  U1AXAT10N  STUDIES... 

KYOHAMBAMAN 


44 


SPECTRAL  HOLE  BURNING  AND  INSTABILITIES  IN 
SEMICONDUCTOR  LASERS 


M.  Sargent  m  and  S.W.  Koch 


PUBUCATIONS  AND  PRESENTATIONS 

W.W.  Qtow,  M.F.  Pereira  and  S.W.  Koch.  ‘A  many-body  treatment  of  the  modulation  response  in  a 
strained  quantum  well  semiconductor  laser  medium,”  Appl.  Phys.  Lett.  61,  7S8  (1992). 

R.  Jin,  D.  Boggavarapu,  G.  Khitrova,  H.M.  Gibbs,  Y.Z.  Hu,  S.W.  Koch  and  N.  Peyghambarian,  *The 
linewidth  broadening  factor  of  a  microcavity  semiconductor  laser,"  Appl.  Phys.  Lett.  61,  1883  (1992). 

M.  Sargent  m,  S.W.  Koch  and  W.W.  Chow,  "Sidemode  gain  in  semiconductor  lasers,"  J.  Opt.  Soc.  Am. 
B  9,  1288  (1992). 

M.  Rose,  M.  Lindberg,  W.W.  Chow,  S.W.  Koch  and  M.  Sargent  m,  "Coroposite-cavity-mode  approach 
to  single-mode  semiconductor  laser  feedback  instabilities."  Phys.  Rev.  A  46,  603  (1992). 

S. W.  Koch  and  W.W.  Chow,  "Microscopic  theory  and  modelling  of  semiconductors  and  semiconductor 
lasers,"  invited  article  for  Trends  in  Optical  Engineering,  publ.  by  Research  Trends,  Council  of  Scientific 
Research  Integration,  India  (1993). 

F.  Jahnke,  S.W.  Koch  and  K.  Henneberger,  "Dynamic  response  of  short-cavity  semiconductor  lasers," 
Appl.  Phys.  Lett.  62,  2313  (1993). 

F.  Jahnke,  K.  Henn^rger,  W.  Schifer  and  S.W.  Koch,  "Transient  nonequilibrium  and  many-body 
effects  in  semiconductor  microcavity  lasers,"  J.  Opt.  Soc.  Am.  B  10,  2394  (1993). 

P.  Ru,  J.V.  Moloney,  R.  Indik,  S.W.  Koch  and  W.  Chow,  "Microscopic  modelling  of  bulk  and  quantum 
well  GaAs-based  semiconductor  lasers,"  Optics  and  Quantum  Electronics  (submitted). 


SCIENTDIC  PERSONNEL 
M.  Sargent  Ed 
F.  Jahnke 


RESEARCH  FINDINGS 

In  this  projea  we  concentrated  on  a  microscopic  analysis  of  semiconductor  gain  media, 
semiconductor  lasers,  and  laser  instabilities  such  as  spectral  and  kinetic  hole  burning.  Given  the  relatively 
small  amount  of  research  funding  through  JSOP,  it  is  probably  fair  to  say  that  this  project  led  to  a  very 
substantial  number  of  results  in  the  semiconductor  laser  arena. 

In  addition  to  the  results  that  were  highlighted  in  our  earlier  research  and  progress  reports,  during 
the  last  part  of  this  project  we  concentrated  on  the  development  and  extension  of  a  truly  quantum 
mechani^  nonequilibrium  theory  of  semiconductor  lasers  with  a  special  emphasis  on  the  novel 
microcavity  laser  structures,  such  as  vertical  cavity  surface  emitting  lasers  (VCSELs)  or  the  microdisk 
whispering  gallery  mode  lasers. 

For  these  microcavity  lasers,  we  investigated  the  influence  of  many-body,  nonequilibrium,  and 
cavity  effects,  and  the  emission  characteristics  and  laser  lineshape.  For  many  operating  conditions  we 
find  that  the  electron-hole  plasma  in  these  devices  is  in  a  genuine  nonequilibrium  state.  The  blocking  of 
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the  carrier  states  in  the  gain  region  leads  to  a  pumping  that  has  maximum  efficiency  in  energetically  high 
lying  states.  Hence,  the  pumping  produces  carriers  with  a  higher  than  average  kinetic  energy. 

At  the  same  time,  the  lasing,  which  takes  place  close  to  the  center  of  the  gain  region,  removes 
carriers  with  less  than  average  kinetic  energy.  The  competition  between  the  different  timescales,  i.e., 
the  carrier  cooling  time  on  the  one  hand  and  the  stimulate  carrier  lifetime  on  the  other,  determines  the 
resulting  nonequilibrium  carrier  distribution.  If  the  stimulated  recombination  time  is  shorter  than  the 
cooling  time,  we  encounter  kinetic  hole  burning  and  substantial  carrier  heating,  which  can  cause 
electron-hole-plasma  temperatures  that  exceed  the  lattice  temperature  by  as  much  as  300  K. 

In  our  investigations  of  the  linewidth  of  the  laser  emission,  we  used  our  self-consistent,  fully 
quantum  mechanical  nonequilibrium  laser  theory  to  compute  the  spectral  emission  characteristics  of 
semiconductor  microcavity  lasers  as  a  function  of  the  pump  intensity.  Our  approach  is  allows  us  to 
predict  the  laser  emission  spectrum  starting  from  the  spontaneous  emission,  well  Mow  threshold,  all  the 
way  to  the  stimulation  dominated  emission,  well  above  threshold.  These  calculations  have  been  done  for 
a  number  of  different  pumping  conditions  (i.e..  injection  pumping  or  optical  pumping  with  various  excess 
energies)  and  different  resonator  conditions.  These  resonator  conditions  influence  the  amount  of 
spontaneous  emission  that  is  coupled  into  the  laser  mode. 

As  one  of  the  most  striking  results,  we  find  (in  agreement  with  experimental  observations  by  R. 
Slusher  et  al.  at  AT&T  Bell  Labs)  that  the  usual  Shawlow-Townes  line  narrowing  for  increasing  laser 
intensity  does  not  occur  in  microcavity  lasers  with  high  spontaneous  emission  coupling  efficiency.  Due 
to  the  substantial  contribution  of  the  spectrally  broad  spontaneous  emission,  for  the  high  coupling  case 
we  find  that  the  laser  linewidth  rapidly  saturates  with  increasing  intensity  and  stays  more  or  less  constant, 
the  spectral  width  being  determined  primarily  by  the  width  of  the  cavity  resonance  under  the  realized 
inversion  conditions. 

In  summary,  our  theoretical  analysis  clearly  demonstrates  the  necessity  to  include  the  quantum 
mechanical  many-l^y  and  Coulomb  effects  into  a  consistent  analysis  of  microcavity  lasers.  Based  on 
this  approach,  a  numM  of  experimentally  relevant  modifications  of  the  laser  characteristics  have  been 
prediaed,  and  (partially)  experimentally  verified.  Even  though  this  project  will  no  longer  be  funded 
under  JSOP,  we  will  continue  our  investigations  to  obtain  a  coherent  understanding  of  the  light-matter 
interactions  taking  place  in  semiconductor  lasers. 
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Pump  Power  (^W) 

Figure  1.  Comperisoo  of  the  experimentel  (open  and  solid  symbob)  and  theoretical  (lines)  results  for  the  linewidth 
and  output  intensity  of  microdisk  lasers.  The  solid  data  poinb  are  for  a  2.2-(aa  Aimmrtrr  microdisk  (^KWlaoeous 
emission  coupling  measured  to  be  0. 1-0.3),  and  the  open  points  are  for  a  S-/un  diameter  microdisk  (spontaneous 
emission  coiqiling  measured  to  be  0.02-0.03).  The  theoretical  results  are  for  spontaneous  emission  couplings  of  0.4, 
0.1,  and  0.01  (from  top  to  bottom). 
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Figure  2.  Effective  carrier  tempenture  versus  pump  rate  for  carrier  injecdoo  pumping  of  a  microcavity  laser  with 
a  cavity  length  of  1  /un  and  effecdve  miiror  reflectivity  of  99%. 
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RESEARCH  FINDINGS 

This  research  was  aimed  at  characterizing  semiconductor  nanostructures  with  an  STM  using 
intense  optical  fields.  We  have  developed  a  novel  approach  of  probing  fast  optical  interactions  using  the 
scanning  tunneling  microscope  (STM).  We  have  demonstrated  the  feasibility  of  this  technique  by  using 
beat  frequencies  of  longitudinal  modes  of  a  HeNe  laser  at  the  tunneling  junction'  by  presenting  the  fast 
optical  response  of  photoexcited  charge  carriers  in  the  layered  structure  semiconductors  n-type  MoS}  and 
p-type  WSej.  These  experiments  open  the  door  to  a  new  and  exciting  field  of  optics  that  takes  advantage 
of  the  short  time  response  of  the  STM  junction,  which  is  several  femtoseconds. 

The  termination  of  our  JSOP  filling  has  put  this  project  in  jeopardy;  however,  we  are  continuing 
to  develop  the  technology  of  locally  probing  the  optical  response  of  carriers  in  quantum  well  structures, 
quantum  dots  and  other  nanostructures.^  Our  laboratory  has  achieved  the  shortest  time  response  of  an 
STM  using  optical  excitations  (2.3  nsec)  of  semiconductors,  and  is  one  of  only  a  few  laboratories  in  the 
United  States  to  perform  such  experiments.  The  experiments  so  far  have  b^n  conducted  in  ambient, 
using  layered  semiconductors.  However,  we  have  just  finished  constructing  a  special  UHV-STM  facility 
that  is  now  fully  operational.  We  are  currently  moving  this  experiment  into  the  UHV-STM  facility.  Our 
specific  plans  call  for  intensifying  die  research  using  the  following  methods:  (1)  incorporating  the 
experiment  into  our  new  optically  integrated  UHV-STM  facility;  (2)  using  an  Ar-pumped  dye  laser  for 
wavelength  scans;  and  (3)  using  the  different  modes  of  the  laser  to  probe  the  temporal  response  of  the 
carriers.  The  samples  to  be  used  will  be  nanolithograpbically  machined  semico^uctor  and 
metal-semiconductor  structures  and  cleaved  multiple  quantum  well  structures. 

A  revised  edition  of  my  book’  Scanning  Force  Microscopy  has  been  written  in  the  course  of  this 
project,  and  will  be  published  by  OxfDrd  University  Press  in  May  1994. 
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RESEARCH  FINDINGS 

Our  recent  research  has  been  in  the  areas  of  femtosecond  pulse  propagation  in  semiconductor 
amplifiers  and  all-optical  switching  using  quantum  dots  semiconductor  materials.  This  work  has  been 
published  in  several  journal  articles,  a  book  chapter,  and  has  been  presented  at  several  scientific  meetings. 

Investigations  of  femtosecond  pulse  propagation  in  semiconductor  amplifiers  has  continued  due 
to  potential  applications  to  all-optical  switching.  In  particular,  the  propagation  of  fs  pulses  is  radically 
different  from  ps  pulses.  For  ps  or  longer  pulses,  gain  saturation  provides  an  adequate  physical  basis  for 
understanding  pulse  propagation:  here  only  amplification  of  die  output  raergy  with  respect  to  the  input 
can  occur.  At  hi^  input  intensities,  the  gain  is  saturated  as  the  carrier  density  approaches  the 
transparency  value  where  gain  equals  loss.  In  the  gain  saturation  model  it  is  assumed  diat  the  electron 
and  hole  distribution  fimctions  remain  Fermi-Dirac  distributions  with  fixed  temperature  but  variable 
density.  This  model  is  therefore  only  valid  if  the  pulse  duration  gready  exceeds  the  carrier  redistribution 
time.  We  have  observed  extreme  dqiartures  from  this  behavior  in  numerical  simulations  in  which  high 
intensity  100  fi  pulses  are  absorbed  in  semiconductor  amplifiers.  This  can  occur  since  the  carrier 
distribution  functions  are  driven  far  from  equilibrium  with  a  large  concomitant  increase  in  the  carrier 
temperature  (from  300  K  to  600  K  in  the  simulations).  This  increase  in  die  carrier  temperature  leads  to 
a  reduction  in  the  gain  and  to  absorption  at  high  input  intensity.  At  a  certain  critical  input  intensity  the 
gain  and  loss  are  balanced  and  the  pulse  propagates  in  a  loss-less  manner,  reminiscent  of  the  solitons  of 
self-induced  tran^arency. 

For  all-optical  switching  iqiplications  a  medium  with  a  large  nonlinear  refractive-index  is  sought 
while  at  the  same  time  minimizing  Ae  detrimental  effects  of  absorption.  In  material  terms  this  calls  for 
a  system  whose  resonance  conditions  can  be  varied  so  that  the  tradeoff  between  refractive-index  and 
absorption  can  be  adjusted  for  a  given  laser  frequency.  Our  recent  work  shows  that  this  degree  of 
flexibility  is  realized  in  quantum  dots  since  the  quannim  confinement  leads  to  a  tunable  blue-shift  in  the 
lowest  exciton  resonance  which  mcreases  with  decreasing  dot  radius.  In  particular,  in  collaboration  with 
Prof.  S.W.  Koch  we  have  calculated  the  all-optical  switching  figure  of  merit  for  a  variety  of  quantum  dots 
materials  using  a  microscopic  semiconductor  model.  These  results  show  explicitly  that  quantum  dots  are 
good  candidates  for  all-optical  switching  iqiplications. 
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RESEARCH  FINDINGS 

Considerable  progress  has  been  made  towards  understanding  the  physics  of  laser  cooled  and 
trapped  atoms.  One  example  is  the  recent  observation  of  quantized  atomic  motion  in  one-,  two-,  and 
thre^imensional  optical  lattices. The  new  and  detailed  understanding  of  the  mechanical  aspects  of 
light/matter  interactions  has  already  lead  to  the  suggestion  of  mechanisms  by  which  atomic  samples  with 
tenq>eratures  in  the  nanoKelvin,  a^  possibly  picoKdvin,  regime  may  be  prepared. 

The  study  of  the  physics  of  ultracold  atoms  opens  up  a  br^  new  field  of  fundamental  and 
applied  research.  It  plays  an  essential  role  in  understanding  die  fundamental  limitations  of  atomic 
cooling,  and  hence  of  devices  such  as  atomic  clocks:  the  many-body  effects  that  are  predicted  to  occur 
when  the  atomic  temperature  is  near  the  recoil  limit  are  bound  to  lead  to  shift  and  broadening 
mechanisms,  and  possibly  nonlinear  heating,  thereby  influencing  the  achievable  performance  of  these 
devices.  Many-body  effects  are  expected  to  be  important  in  lithographic  applications  as  well,  whose  basic 
principles  are  closely  related  to  those  of  optical  lattices:  to  achieve  reason^le  writing  speeds,  high  atomic 
densities  will  clearly  be  required. 

In  addition  to  these  immediate  practical  concerns,  the  availability  of  ultracold  atoms  is  of 
considerable  interest  in  atom  interferometry,  since  it  naturally  leads  to  die  possibility  of  coherendy 
splitting  atomic  wave  functions  with  a  large  s^aruion  between  the  arms  of  the  interferometer.  Such 
macroscopic  separations  are  highly  desirable  since  most  applications  of  atom  interferometry  require  that 
the  arms  of  die  interferometer  be  placed  in  different  environments  and/or  diat  large  areas  be  covered  by 
the  interferometers,  as  e.g.  in  gyroscopic  applications.  Cold  atoms  also  have  profound  implications  in 
the  physics  of  collisions,  since  they  lead  us  to  a  regime  where  standard  approximations  cease  to  be  valid.’ 

The  net  effect  of  any  two-body  interaction  V^2(R)  is  given  by  its  overlap  with  the  atomic  wave 
functions,  properiy  symmetrized  in  the  case  of  indistinguishable  particles.  For  ^ort-range  interactions, 
a  sizeable  overly  bc^een  the  atomic  wave  functions,  and  hence  a  high  atomic  density,  is  necessary  to 
achieve  noticeable  effects,  but  in  the  case  of  the  long-range  dipole-dipole  interaction  b^een  coherendy 
driven  atoms,  many-body  effects  can  be  expected  already  at  relatively  modest  densities.  In  this  situation, 
it  becomes  incorrect  to  treat  the  atoms  in  die  sample  as  independent.  Rather,  many-body  effects  lead  to 
an  effective  nonlinear  bdiavior  of  die  single  atoms.  This  regime,  which  we  call  'nonlinear  atom  optics,' 
is  to  de  Broglie  optics  what  nonlinear  optics  is  to  conventional  optics.  It  suggests  that  novel  effects  such 
as  the  generation  of  atomic  solitons  and  solitary  waves,  atomic  wave  mixing,  and  atomic  phase 
conjugation  could  be  possible.  In  addition,  sufliciendy  high  atomic  densities  in  optical  lattices  could  lead 
to  spontaneous  pattern  formation  and  the  creation  of  other  stable  long-range  structures,  such  as  possibly 
superlattices.  Even  more  intriguing,  perhaps,  is  that  the  combination  of  nonlinearities  with  dissipation, 
in  the  form  of  ^ntaneous  emission,  and  of  a  'pump'  mechanism,  might  suggest  the  possibility  of  a 
'coherent  atomic  beam  generator,'  which  could  loosdy  be  considered  as  the  atom  optics  equivalent  of 
a  laser. 

During  this  year,  we  have  initiated  a  theoretical  program  to  study  diese  effects.  This  is  a  high- 
risk,  high-pay-off  task  in  nearly  virgin  tenitory,  where  two  major  avenues  of  investigation  are  possible. 
The  first  one  consists  of  concentrating  one’s  efforts  on  the  understanding  of  the  dynamics  of  the  system 
at  the  microscopic  level.  In  practice,  two  atoms  are  about  as  much  as  can  be  handled  in  this  approach, 
which  reduces  then  to  the  theory  of  diatomic  collisions  between  ultracold  atoms.  This  challenge  is 
presendy  being  tackled  in  particular  at  NIST,  in  the  group  of  Dr.  Julienne. 

We  have  chosen,  instead,  a  more  'phenomenological*  approach,  paying  relatively  little  attention 
to  the  detailed  way  the  atoms  interact  at  the  microscopic  level,  and  concentrating  on  more  'global*  effects 
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such  as  the  generation  of  solitons  and  spatial  patterns,  phase  transitions,  etc.  In  the  absence  of  a  deuiled 
understanding  of  cold  atom  collisions,  such  an  approach  has  considerable  merit  in  permitting  one  to 
investigate  the  possible  effects  of  these  two-body  interactions  in  a  sample  of  ultracold  atoms.  In  this  case, 
the  specific  form  of  the  two-body  interaction  is  considered,  in  a  sense,  as  a  ’ftee  param^er,'  and 
predictions  about  the  behavior  of  the  system  can  be  made  for  various  forms  of  the  potential. 

The  general  many-body  Hamiltonian  including  a  two-body  interaction  can  be  expressed  in  a 
standard  way  as 


♦  1 1  dld2d3d4{\ ,2 1  V\  2A)^\\)PO)k^)ky) 

-  ♦  V  . 

Here,  is  the  single-particle  Schrddinger  Hamiltonian,  and  includes  the  atom-field  dipole  interaction, 
the  electromagnetic  field  being  either  classical  or  quantized,  and  Vis  the  two-body  dipole^ipole  potential. 
All  kets  1 1)  give  a  complete  description  of  the  single  atom  states,  i.e.  they  include  both  the  center-of- 
mass  and  the  internal  quantum  numbers.  We  use  numbers  to  indicate  dummy  variables  and  letters  to 
label  a  specific  state.  The  state  |  f )  is  obtained  bv  application  of  the  creation  operator  kif)  to  the  atomic 
vacuum  |0)  ,  ^^(f)|0>  ■  |/)  ,  where  ■  S(l-l')  for  bosons,  and  the  5-fonaion  should 

be  interpreted  as  a  product  of  ^-functions  for  continuous  quantum  numbers  and  Kronecker  deltas  for 
discrete  ones.  A  normalized  Af-particle  state  1^)^  is  obtain^  as 

I ^  -  -tL  f  . N)km...k(l)IO>  ,  (2) 

where  the  N-body  wave  function /y(l,2,...,A0  is  totally  symmetric  in  its  arguments. 

In  gener^,  the  many-body  problem  described  by  the  Hamiltonian  (1)  cannot  be  solved  exactly. 
There  are,  however,  a  number  of  aq)proximate  techniques  that  have  been  developed.  For  the  bosonic 
problem  at  hand,  time-dq)eQdent  Hartree  techniques  originally  utilized  in  nuclear  physics  seem 
particularly  attractive.  The  central  tenet  of  the  Hartree  approximation  consists  in  factorizing  the  A^-body 
wave  function /^l,..., AO  as  a  product  of  effective  single-particle  wave  functions.  In  the  time-dependent 
Hartree  approximation,  foe  equations  of  motion  for  foe  effective  single-particle  states  are  determined  from 
foe  Hartree  variational  principle.*  This  yields  foe  system  of  effeaive  single-particle  nonlinear  Schrddinger 
equations 

(3) 

♦  (W-l)  pMMJ  <f  .1 1  V|2.3> 

Equations  (3)  form  foe  basis  of  nonlinear  atom  optics.  They  show  that,  from  the  perspective  of  a  single 
atom,  the  N  -  I  other  bosonic  atom  effectively  aa  as  a  nonlinear  medium.  This  is  analogous  to  the 
situation  in  conventional  nonlinear  optics,  where  foe  presence  of  a  medium  can  lead  to  an  effective 
nonlinear  b^avior  of  foe  light  field  when  foe  medium  dynamics  is  traced  over. 

Equations  (3)  are  quite  general,  as  the  single-particle  Hamiltonian  Ho  could  describe  the 
interaction  between  atoms  a^  a  standing-wave  field,  their  interaction  with  a  magneto-optical  trap,  the 
potential  of  an  optical  lattice,  or  any  situation  that  does  not  involve  many-body  effects.  Also,  foe  atoms 
could  be  simple  two-level  or  multi-level  systems  as  are  required  in  most  standard  laser  cooling  schemes. 

We  have  already  analyzed^  a  number  of  consequences  of  Eqs.  (3).  In  particular,  for  foe  case  of 
two-level  atoms  interacting  with  a  near-resonant  light  field,  the  so-called  near-resonant  Kapitza-Dirac 
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effect,  we  predict  modtficatioiis  of  the  conventional  PendellOsung  of  Bragg  scattering.  In  another  regime, 
the  problem  is  reduced  to  a  classical  massive  Thirring  model,  with  the  possibility  of  generating  atomic 
Thirring  solitons.  In  one  dimension,  gap  solitons  as  well  as  nonlinear  focusing  effects  are  expected  to 
occur.  Furthermore,  if  the  equations  of  nonlinear  atom  optics  are  transformed  to  momentum  space  they 
reveal  a  variety  of  four-wave-mixing  terms  including  those  responsible  for  phase  conjugation.  This 
suggests  the  possibility  of  an  atomic  phase  conjugator  capable  of  reconstructing,  or  time-reversing,  atomic 
wave  functions. 

Turning  now  to  the  question  of  numbers,  we  note  that  nonlinear  effects  are  expected  to  become 
noticeable  as  soon  as  the  effective  nonlinearity  (^-l)Vg|^|Vq  becomes  comparable  to  the  effects  of 
linear  diffraction,  whose  strength  is  of  the  order  Stated  differently,  this  indicates  that  nonlinear 

effects  should  arise  as  soon  as  the  two-body  interaction  energy  becomes  comparable  to  the  photon  recoil 
energy.  For  the  long-range  dipole-dipole  interaction  and  at  resonance,  this  yields  the  critical  linear 
density  a>  NIL  a  w^ql0y.  Here  /3  is  a  geometrical  factor  that  reflects  the  fact  that  the  effective 
strength  of  the  dipole-dipole  interaction  is  given  by  the  fraction  of  spontaneously  emitted  photons  that  can 
be  reabsorbed  by  atoms  inside  the  laser  beam.  For  the  Cesium  Dj-line,  we  have  p^  a  30//3  cm'‘,  which 
corresponds  to  relatively  modest  densities  of  p  at  3  *  10'°  cm  **  for  -  10*’.  Clearly,  these  numbers 
are  consistent  with  the  experimental  observation  that  cooling  ceases  to  follow  single  atom  scaling  laws 
for  the  densities  diat  can  now  be  achieved  in  the  lab.*'*''°  A  very  important  point  here  is  that  nonlinear 
effects  are  expected  to  occur  much  before  the  occurrence  of  Bose  condensation,  and  hence  experiments 
are  possible  without  major  breakthrou^  in  cooling  and  trapping  techniques. 

The  time-dq)endent  Hartree  theory  considers  only  the  coherent  part  of  the  many-body  problem. 
Of  course,  as  is  the  case  in  conventional  atom  optics,  the  incoherent  effects  of  spontaneous  emission  are 
expected  to  destroy  the  collective  effects  discussed  above  to  some  extent.  We  have  now  staned  to  study 
the  conq)etition  b^een  coherent  and  incoherent  effects  as  soon  as  we  have  a  solid  understanding  of  the 
coherent  effects  alone.  This  is  an  essential  step  both  from  a  practical  and  a  fundamental  point  of  view. 
In  particular,  atomic  solitons  are  mesoscopic  quantum  systems,  which  will  permit  us  to  study  the 
transition  from  the  microscopic,  quantum  world  to  the  macroscopic  classical  world.  It  is  now  well 
understood  that  the  survival  of  such  systems  ostt  a  reasonable  time  scale  depends  critically  on  their 
coupling  to  the  environment.  In  the  present  case,  this  coupling  is  mainly  via  spontaneous  emission. 

The  description  of  die  combined  coherent  and  incoherent  effects  is  best  performed  in  terms  of 
a  master  equation,  which  includes  both  die  effects  of  the  vacuum  field  and  of  the  laser  fields  with  whic 
the  atoms  are  interacting.  Specifically,  one  needs  to  derive  an  equation  of  motion  for  the  effective  single 
particle  density  operator  p(2,l)  ■  s^^(l)^(2)>  •  Tr[^\l)^{2)p],  where  p  is  the  Afiiarticle  density 
operator.  This  equation  can  be  derived  by  applying  the  Bom-Markov  approximation  for  the  vacuum 
modes  of  the  reservoir,  and  the  Haitree-Fock  approximation  to  factorize  higher-order  correlation 
functions  of  the  form  into  products  of  second-order  correlation  functions,  thereby 

introducing  direct  and  exchange  contributions.  The  resulting  master  equation  contains  the  coherent 
contributions  already  present  in  the  pure  state  Hartree  approach.  In  addition,  it  contains  a  dissipative  part 
that  accounts  for  spontaneous  emission  and  the  concomitant  atomic  recoil  as  well  as  a  nonlinear  non- 
Hamiltonian  term. 

To  properly  assess  the  impact  of  spontaneous  emission  on  the  many-body  effects  that  we  discussed 
earlier,  it  is  necessary  to  carry  out  numerical  simulations,  which  are  now  underway.  We  also  have 
initiated  an  effort  to  check  die  results  of  the  effective  single  atom  Hartree-Fock  approa^  against  "exact* 
results.  For  this  reason,  we  also  carry  out  explicit  calculations  with  a  small  number  of  atoms,  say  two 
or  three,  using  Monte  Carlo  wave  function  simulations. 
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